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We   have   updated  our   analytic  representations   of 
fundamental  molecular  cross  sections  for  H2C  vapor,   drawinq 
extensively  on  the  recent  work   of  the  Applied  Physic  Group 
at   Los  Alamos.     These  cross   sections,   with   slight 
modifications,  have  been  used  in  a  Monte  Carlo  computer  code 
to  build  up  a  statistical   profile  of  the  collisional  events 
incurred  by  an   electron  as  it  is  deqraded   in  H2O.    The 
results  are  described   by  a  continuous  function   called  the 
spatial  yield  spectrum.    This  four  dimensional  function  can 
then  be  used   to  calculate  the  spatial   distribution  of  the 
excited  and  ionized  states  created  by  the  electrons  which  in 
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turn  are   used  to  describe  tke   steady  state  input   into  a 
Diffusion-Kinetic  computer  code.     In  addition,   reaction 
rates   and   diffusion   coefficients   which   are   not 
experimentally  known   have  been   estimated  usinq   current 
models.    It   is  then  possible   to  calculate   the  spatial 
distribution  of  the  final  products   for  a  steady  state  beaa 
of  electrons.     Additional  results  related  to   the  yield 
spectruia  such  as   eV/ion  pair  and  deposition   contours  are 
also  presented-    The  results  are  intended  to  facilitate  an 
understandinq  of   the  effects   of  radiation   in  bioloqical 
systems. 


IX 


CHAPTER  I 
INTRODUCTION 


Overview 
Radiation  physics,    radiation  chemistry   and  radiation 
biology  form  a  triad  of   scientific  fields  which  address  the 
effects   of  radiation   in   terms   of  three   basic   and 
historically  different   points  of   view.    The   stady  of 
radiation  effects   on  bioloqical  systems   has  in   the  past 
looked  toward  chemistry  to   supply  the  fundamental  concepts 
necessary  to  understand  this  field.   Radiation  chemistry  nas 
been  able   to  explain   many  of   the  observed   effects  of 
radiation  with  the  free  radical   theory  and  the  discovery  of 
the  aqueous   electron.    At  the   same  time  the   field  of 
radiation  physics  has  evolved  from   the  study  of  macroscopic 
concepts  suca  as   exposure,   dose  and  LET   (linear  enerqy 
transfer)   to  the  study  of   individual  atomic  and  laolecular 
processes  which  require  the   formalism  of  quantum  mechanics. 
Durinq  the  past  few  years  the   proqress  in  these  fields  has 
increased  to   the  point  where  it   is  possible  to   beqin  to 
understand  the  bioloqical  effects  of   radiation  in  terms  of 
individual  atomic  events  and  their  aqqreqate  behavior!   This 
is  due   in  large   part  to  the   development  of   high  speed 
computers,    experimental  studies   of   basic  atonic   and 
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molecular  cross  sections  and   picosscond  studies  of  chemical 
kinotics. 

While  the  incident  radiation  may  be  qamma  rays,   protons 
or  charged   particles,   ultimately  it  is   electrons  which 
dissipate  much  of  the  incident  enerqy.   Radiation  deposition 
is  split  into  two  broad  categories:   (1)  direct  effects--  in 
which  the   incident  radiation  interacts  directly   with  the 
bioloqically  sensitive  sites  and  (2)   indirect  effects--  in 
which   the  incidejit   radiation  interacts   with  the   laaior 
component  of  the  system,  for  example  H  0  which  in  turn  forms 
free  radicals  which  interact  with  the  important  biocheaicals 
(see  fiq.  1-  1)  . 

Referinq  to  this   figure  let  us  describe   the  knowledge 
necessary   to   understand  the   fundamental   physics   and 
chemistry  which  leads  to  biological   damage  due  to  incident 
electrons.    First,   detailed  knowledge  of  the  medium  into 
which  the  electrons  will  be  fired  is  required.    This  would 
imply  knowledge  of  the  various   molecules  which  populate  the 
medium,  their  spatial  relationship  and  individual  atomic  and 
molecular  energy  levels.    Next,  ccllisional  cross  sections 
for  excitation,   ionization  and  dissociation  for  individual 
molecules   are  needed.     These  include   both  total   and 
scattering  cross  sections.    At  this   time  we  have  spatial 
distributions  of  "initial  yields"  of,  for  example  H  C+,  0H+, 
H+,  OH,   H  and  H  0*  (excited  states).    Part  oi  the  radiate 
enerqy  has  oi  course  been   absorbed  directly  by  the  critical 
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biochemicals.    We  move  now  from   the  physical  statje  to  the 
chemical  physical  staqe.   Furthermore,  since  the  bulk  of  the 
dauiaqe  is  caused  by  indirect  effects  additional  chamical  and 
physical  properties  are  required.   A  larqe  percentaqe  of  the 
information  needed  in  this  study  is  drawn  from  literature  in 
the  combustion  field  (Kee  et  al- ,  19  83:   Sestbrook,   198  1) 
since  hydroqen-oxyqen  flames  form  a   maior  area  of  study  in 
this  field  and  many  of  the   chemical  species  present  in  this 
type  of  flame  are  also  present  in  the  irradiated  H,^0  system. 
The  ions  and  excited  states   quickly  react  or  dissociate  and 
leave  us  with  a  distribution  of  neutral  products.   Therefore 
while  it  is   necessary  to  know  the   possible  reactions  and 
associated  reactions  rates  for  the  ion  and  neutral  products, 
knowlodqe  of  the  diffusion  coefficients   for  the  ions  is  not 
as  critical  as   that  for  the  neutrals.     Lastly,   these 
molecules  are  free  to  move  (within  bounds)   and  therefore  it 
is   necessary  to   know   transport   properties  for   them. 
Ultimately  it  is   necessary  to  know  the   types  of  chemical 
reactions   which  result   from  the   interaction  of   these 
products  and  the  nearby  biomoleculos.    This  information  is 
found  in  the  radiation   chemistry  literature  (Scholes,   et 
al. ,  1969;  i-'yers,  1974)  but  often  does  not  take  into  account 
bioloqicai  structure  and  function  which   will  of  course  lead 
to  modifications   of  the  immediate  chemical   effect.    For 
exaffi:ple,   it  is  well  known  that   DNA  will  repair  itself  if 
damaqed.   Furthermore  the  ability  to  repair  itself  is  also  a 
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function  of  the   amcunt  and  or  type  of   damaqe  sustained  by 
the  molecule. 

While  it  is  not   possible  at   this  time   to  study  a 
realistic  biological  system  iu  this  manner,   several  studies 
usinq  this  method   have  been  applied  in   varyinq  daqrees  to 
study  the  radiolysis  of  H^O  in   both  its  liquid  and  qaseous 
state.   The  ma-jor  proqrams  are  located  at  Oak  Ridqe  National 
Labs  (Hamm  et  al.  ,  1978;   Sitchie  et  al. ,  1978),  Los  Alamos 
National  Labs  (Zaider,   Brenner   and  Wilson,   1982  (ZBW) ; 
Zaider  and  Brenner,   1983;    Brenner  and  Zaider,   1983), 
Lawrence  Berkeley  Laboratory  (Maqee   and  Chatteriee,   1980; 
Chatterjee  and   liaqee,   1930)    and  Battelle   Northwest 
Laboratories  (Wilson  et  al. ,   1978).    Much  of  the  work  by 
these  qroups  depends  heavily  on   the  pioneerinq  work  done  on 
fundamental  collisioaal  cross  sections  of  electrons  with  H 2O 
at  the  University  of  Florida   (OF)   (Green  et  al. ,   1971; 
Olivero   et  al. ,    1972;   Kutcher   and  Green,    1976). 
Subsequently  the  thrust  of  the   research  at  U7  has  involved 
usinq  similar  cross   section  work  to  study   the  collisional 
history  of  electrons  in  various   qases  and  representinq  this 
information  in  analytic  form   (Green,   Garvey  and  Jackman, 
1977;  Jackman,  1978;  Green  and  Siuqhal,  1979  (GS) ;   Sinqhal 
and  Green,   1981;   Sinqhal  et  al. ,   1983).    Kany  of  the 
techniques  required  for  these  calculations  oriqinated  in  the 
field  of  atmospheric  physics.    With   this  thesis  we  brinq 
toqether   the  work   on   collisional   cross  sections   and 
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transport  calculations  at  UF  to   iaply  the  initial  yields  of 
ionized,   excited  and  dissociated  H  2O  in  tiie  qaseous  state. 
These  are  then  used  as   input,   alonq  with  chemical  kinetic 
rate  constants  and  diffusion   coefficients  to  calculate  the 
spatial   distribution  of   chemical  species   created  by   a 
iaonoenerqetic  beam  of  electrons. 

Hi s to r Y  and , Mo tiyati on 
The  field  of  radiation  chemistry   originated  in    the  18th 
century  with   the  study  of  chemical  cJianqes   produced  by 
electrical  discharqes   in  qases.     Many  of   the  reaction 
mechanisms  to   explain  the  experimental  results   studied  in 
the  iaterveninq  years  still  remain  opaque  (Anderson,   1968) . 
The  establishment  of  radiation  chemistry  as  an  exact  science 
is  basically  due  to  the  work  of  Lind  (192  1,   1930)   in  the 
U.S.A.   and  Muud  (1935)  in  Belqium.    Much  01  this  research 
used  alpha  particles  as  a  radiation   source  and  many  of  the 
quantitative  measurements   are  still   qood  today.     For 
example,   it  vjas  observed  that  the  radiolysis  of  pure  water 
vapor  produced   almost  no  siqnificant   decomposition  unless 
chemically  reactive   scavengers  were   introduced  'into   the 
system.     The    initial   interpretation   of  the   observed 
chemistry   was   centered   on  the   assumption   that   the 
r adioloqically  produced  ions  were   surrounded  by  "clusters" 
of  neutral  molecules  which  would  ultiiaately  recoiubiae  either 
with  themselves,   if  they  were  of  opposite  charqe,   or  with 
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free  electrons,    Dpon  doinq  so  the  enerqy  released  would  be 
shared   by  the   adjacent  neutral   molecules  which   would 
instigate  further  reactions. 

as  early  as   1909  Debierne  (1909)   proposed   that  free 
radicals   were  responsible   for   the  observed   chemistry. 
During  the  1930 's   two  papers  by  Eyrinq,    Hirshf elder  and 
Taylor  (1936a,   1936b)   indicated  that  since  the  "W  value" 

(average  energy  to  form  an  ion   pair)   is  much  greater  than 
the  lowest  ionization  energy  of   a  laolecuie  that  the  excess 
energy  must   be  channeled  directly  into   excitation.    In 
addition  thay  proposed  that   charge  neutralization  reactions 
would  also  lead  to  excited  states.   The  excited  states  would 
then   dissociate  forming   free  radicals   from  which   the 
remainder  of  the  chemistry  would  follow. 

Currently  it   is  realized  tnat   both  ions   and  excited 
molecules  produced   directly  by  the  incident   radiation  are 
important,    Recent  research  by  Anderson,  Knight  and  Winter 

(1964,1966a,  1966b)  has  clarified  the  neutralization  process 

in  H  O.     These  reactions   produce  free   radicals  which 

2 
continue   to  build   up   until   their  concentrations   are 

sufficiently  high   tnat  reverse  cheaiical   reactions  compete 

with   radical   production   and    small   steady   state 

concentrations  of  H  ,  O  ,  and  ii  0   are  observed. 

2    2        2  2 

The  choice  of   using  fl  0  vapor  instead  of   liguid  as  the 

2 
first  test  case  for  study   reguires  a  coaraent.    While  the 

radiolysis  of  H  0   liguid  is  obviously  of   more  biological 
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significance  tlian   H  2O  vapor,   physical  processes   in  the 
liquid   are  not  well  known.    For  example,   the  question  of 
defininq  ionization  in  condensed  matter  (Douthat  1983)   does 
not  need  to  be  addressed.    In  coaparison,   the  collisional 
cross  sections  for  the  qas  phase  are  known  to  a  Quch  qreater 
deqree  of  accuracy.   Often  many  of  the  liquid  cross  sections 
are  extrapolated   from  qaseous  models  (Kutcher   and  Green, 
1976;   Hamin  et  al.  ,   1978,   1983)   and  questions  of  their 
validity  have  been  raised  by  luoXuti  (1983).   Furthermore  as 
expected  the  general  features  of  both  vapor  and  liquid  water 
are  very  similar  (Thomas,   1969).     The  initial  yields  of 
ionized,  excited  and  dissociated  620  are  comparable  (Thoicas, 
1969;   Kutcher  and  Green,    1976).    However  the  reaction 
mechanism  for  the  liquid  is  more  complex  due  to  the  presence 
of  the   hydrated  electron,    theoretically  predicted   by 
Platzman  (1955)  and  experimentally  observed  by  Boaq  and  Hart 
(1963)   and  the   more  ambiguous  role  of   excited  molecules 
(Thomas,   1969).    Finally  H2 C  vapor  represents  a  system  in 
wnich  the  sensitivity  to  cross  section  inputs  is  most  easily 
studied,   comparisons  witii  experimental  data  are  inherently 
simpler  and   the  calculations   accomplished  for  'the  vapor 
phase   can  be   modified  for   tae  liquid   phase  once   the 
appropriate  input  data  are  obtained. 
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£jgscri£jtioii 
The  purposes  of  this  dissertation  are  to 

1.  Describe  multiple   elastic  scattering  of   low  enerqy 
electrons  in  H^O  vapor   and  characterize  the  spatial 
distribution   of   the  resultant   distribution   of 
electrons. 

2.  Extend  the  work   on  initial  yields  due   to  incident 
electrons  in  H_0  vapor  as  reflected  by  the  concept  of 
a  yield  spectrum  (see  caapter  II) 

3.  Use  the   above  information  to  study   the  subsequent 
chemical  reactions  which  occur   durinq  the  diffusion 
of  these  chemical  species  in  H 2O  vapor. 

This  will  be  developed  in  three  maior  seqments.    In  chapter 
II,   using  the   earlier  work  by  Kutcher   and  Green  (1976) 
multiple   elastic  scattering   distributions  for   electrons 
undergoing   only  elastic   collisions   in   H _0  vapor  are 
calculated.   These  results  along  with  a  compilation  of  cross 
sections  based  on  empirical  models  developed  by  Green  et  al. 
(1971),  Olivero  et  al.  (1972)  and  Zaider  et  al.   (1982)  are 
used  to  calculate  spatial  yield   spectra  (Green  and  Sinqhal, 
1979).    These  are  then  used  to  calculate  spatial  yields  for 
all   the  important   "initial   products"   and  a   spatial 
distribution   for   subexcit ation   electrons.     Analytic 
representations  of   all  these  quantities  are   given-    In 
chapter  III  the  above  calculated   quantities  are  followed  as 
they  diffuse  and  react.    This  is  accomplished  by  solving  a 
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system  of  nonlinear  partial   differential  equations.    This 
approach  is   based  on  the   "classical"  diffusion   model  as 
developed  for  the   radiolysis  of  H^O  which   oriqinated  with 
the  radical   theory  of  Weiss   (1944)   and   expandinq  spur 
kinetics  as  hypothesized   by  Lea  (1946,   1947)    and  was 
ultimately  formulated   by  Safisuel  and  Maqee   (1953),    The 
computations  were   developed  by  Ganquly  and   Maqee  (1956), 
Flanders  and  Fricke  (1958),   Dyne   and  Kennedy  (1958)   and 
Kuppermann  (1961  ,1974).    Recently  these  calculations  were 
revived  by  Magee  and  Chatteriee   (1930)   and  Chatteriee  and 
Maqee  (1980),    Essentially,  it  consists  of  representinq  the 
"initial"  (in  this  instance  those  chemical  species  which  are 
seen  approximatly   10  (-13)   seconds   after  the   radiation 
impact)  yields  of  radicals  by  a  few  estimated  parameters  and 
followinq   their  evolution   in   time,   assuming   spatial 
variation  in  only  one  dimension,   radial  in  either  spherical 
or  cylindrical  coordinates,    depending  on  the  L2T   of  the 
radiation.    The   calculations  are  then   compared  against 
experimental  G-values(#  pcoduced/100  eV)   of,   for  example. 

In  this   study  a  maior  step   is  taken  to   elucidate  the 
mechanisms  involved   in  this  transition  reqion   between  the 
physical  and  early  cnemical   physical  stages.    The  problem 
will   be  solved   in   cylindrical  coordinates,    assuming 
cylindrical  symmetry  and  in  the  steady  state  limit.   Solving 
the  temporal  problem   usiuq  this  technique  would   have  been 
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excessively   expensive.    The   maior   emphasis  of   this 
cofflputatioa  is  to  establish  initial   spatial  yields  for  the 
chemical  species  produced  by  the   incident  radiation  and  to 
tie  the   physics  and   chemistry  toqether   by  taStinq   into 
account  the  ion  chemistry  as   well  as  radical  recomaination. 
The  calculations  are  performed  usinq  the  IMSL  finite  element 
partial  differential  equation  solver  TWODEPEP  (1981).    The 
systeia  of  equations  is  solved  in   the  steady  state  case  and 
spatial   concentrations   of  all   species  involved   are 
calculated.    Finally,    in  chapter   IV  the   quantities 
calculation  are  appraised  and  applications  and  extensions  of 
the  computations  are  discussed. 


CHAPTER    II 
THE    PRODUCTION    Of    INITIAL    YIELDS 


Tae_Mgn t e_Ca r lo„ Me th  gd 

The   Honte  Carlo   approach  to   solvinq  the   transport 
equation  is  a  stochastic  method  of  imitating  the  actual  path 
of   an  electron   as  it   is   deqraded  in   enerqy  by   its 
interactions  with  the  mediutn,     Durinq  the  siaulation  the 
moveaent  of  the  electron  is   governed  by  its  cross  sections 
which  describe   elastic  or   inelastic  collisions   and  the 
generation  of  a  few  random  numbers.   These  numbers  take  into 
account  the  statistical  nature   of  the  degradation  process. 
Depending  on  the  incident  energy   of  the  electron  the  number 
of  incident  electrons  allowed  to   degrade  is  varied  in  order 
to  build  up  a  statistical  profile   of  the  actual  solution  of 
the  transport  equation.  The  statistical  error  resulting  from 
the  Monte  Carlo   calculation  is  discussed  in   the  thesis  by 
Jackiaan  (1973)   where  a  detailed   description  of  an  earlier 
version  of   the  program  used   in  this  calculation   may  be 
found.    A  compact  overview  of  this  particular  study  may  be 
found  in  Jackman  and  Green  (1979), 

In  brief,   the  information  necessary   to  carry  out  this 
calculation  is  a  detailed  collection   of  atomic  or  molecular 
cross  sections,  the  attributes  of  the  incident  electrons,  in 
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particular  their  energy  and  direction  and  the  density  of  the 
mediuia  which  is   arbitrary  since  the  collision   pattern  is 
simply  scaled  spatially  by  this  quantity.   Thus  the  ranqe  of 
the  electrons  which  depends  on  the  density  is  chosen.    The 
information  generated  at  each  collision  for  each  electron  is 
stored  on  magnetic   tape  allowing  data  analysis   to  be  more 
flexible,   since  new  items  of  interest   nay  be  pursued  at  a 
later  time  without  rerunning  the  program.    The  information 
stored  on  the  tape,  at  each  collisional  point,   is  the  type 
of  collision   and  associated   energy  loss,    the  absolute 
position  of   the  electron  collision   and  the   incident  and 
final  energy  of  the  electron,   its  absolute  direction  and  a 
tag  delineating  whether  it  is   a  primary  or  later  generation 
electron.    Extensive   Eodif ica tions  of   the  Monte   Carlo 
proqraa  have  been  jaade  while  applying  it  to  K  O.    Those  are 
documented  in  appendix  A. 

Basic_Molecu,lar_Cross__Sectioas 
lis stic.jCross_ Sect  ions 

The  production   of  a   spatial  yield   spectrum  requires 
detailed  comprehensive   data  on  the  elastic   scattering  of 

electrons  from  H  0.   Unfortunately,  such  data  are  sparse  or 

2 

nonexistent.   To   overcome  this  problem   we  rely   on  the 

extensive  work  done  by  Porter   and  Jump  (1978),   which  was 

recently  updated  at  Los  Alamos   National  Labs  (ZBN)   (1962) 

for  H  0,    to  characterize   elastic  cross   sections  with 
2 

convenient  analytic  functions. 


14 

For  electron   scattering  in  ii  ^C    we   use  the   standard 
Rutherford  representation   for  differential   elastic  cross 
sections  modified   by  a   screeninq  parameter   for  electron 
energies  above  200  eV  (Jackman,  1978).    The    functional  form 
is 

d0(e.E)  =  Z^e^ 2.1 

m  V  (i+2n-cos  e) ^ 


where  r,  (E,Z),tiie  screeninq   parameter,    is  a   modified 
formulation  proposed  by  Grosswendt   and  Waibel  (1973)   and 
used  by  Zaider  et  al.  (1982)  for  H2O.   It  is  qiven  oy 

n   =  n^k/[E(E/mc^+  2)] 

n^  =  a  +  plnE  2.2 

2/3 
k  =  kj^mz"^ 


where  Z  is   the  atomic  number  of  the  medium   and  e,p,v  and 
theta  are   the  charge,   moiaentua,   velocity   and  angular 
deflection  of  the  electron.    The   parameters  are  listed  in 
TaDle  2.  1. 

Below  200  eV   the  empirical  function  proposed  'by  Porter 
and  Jump  (1978) 


dc(O.E)oc  1 +  P(E) , 

dn        1  +  2\(E)  -  cos  e    [1  +  26(E)  +  cos  6]" 


2.3 
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TABLE  2.  1 
Differential  and  total  elastic  scattering  parameters. 


Elastic  Scattering  (  <  200  eV  ) 
PAX 


7.51525        2.9612        -1.7013 

•0.419122      -0.26375       -1.48284 

0.72017     (-2)       0.4307       (-2)        0.6331 

•0.4645       (-4)     -0.26895     (-4)     -0.10911 

-0.10289     (-6)       0.583505(-7)        0.8358       (-2) 

-0.2388       (-3) 

-3.32517 

0.10996 

-0.45255    1-2) 

0.58372     (-4) 

0.58372    (-4) 

-0-24659     (-6) 

0.24775    (-1) 

-0.296264  (-4) 

-0.  120655(-b) 

Eutlierford    Scattering     (    >    200   eV    ) 


=    1.64  P        =    .0825     ^x       =    .000017 

Total   Inelastic    Cross   Section 


Eq  =    1.798  Fj^   =    72.53         ?2   =    11.093 

X      =    -.77  Gi    =    .4734         G2    =    19.9223 


U      =    .00195         Ej^    -    .172  E2    =    13.93 

V      =    150. 
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is   used    wiiece  p  ,    5      and     x     ^^^   parameters    fitted   by   Brenner 
(1982)       to   experimental   data.         These  functional  forms  are 
given   in    Table    2.  2, 

The    experimental    data    used    in   constructinq    this   model    are 
those    oy    Seng     (1975)     (0.35-10    eV,    20-108    degrees),      Traiiaar 
et    al.        (1973)     (15-53    eV,     10-90    deqrees)  ,       Nishimura    (1979) 
(30-200   eV,     10-140    degrees),    Hilqaer    et    al-       (1969)     (60-300 
eV,    35-150   degrees),    and    Bromberq    (1975)     (300-500    eV,       2-160 
degrees) .      Examples   of    the   functions  along    with   experimental 
data   are    presented    in    Figure   2.1, 

For   the      total    elastic      cross    section      we    again      use    a 
function    suggested    by    Porter   and    Jump    (1978)       and    fitted    to 
experimental    data    by    ZBK    (1982) 

„x  2  F   gI 

n(n+l)[V        +E        ]      ^^^      ^^-^n^    "^     n 

where    n=U/3      and   the   parameters     are    listed    in      Table    2.1. 
Data    used   in       fitting    this    function  come      from   Senq     (1975), 
Nishimura     (1979),         Bromberq    (1975)         and    Bruche       (1929) 
(U-36eV).         A    plot      of    the    total    elastic      cross   section   is 
presented    in    Figure    2.2w 
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TABLE  2.2 
Differential  elastic  scattering  functions. 


ill  Pi-1^'"'' 
P(E)=e 

6(E)=e 

«        i-l 

ill'^i-l^ 
X(E)=:e  .35eV<E<10eV 


lo       i  -fi 
.^   \.E  ^  ^ 

=e   ^^   ^  10eV<E<100eV 


13     i -1  1 

illl^i^ 
=  e   ^  -^J-  ^  100eV<E<200eV 


EICRGY-0.60 


-g^  5J3m  mTm  9o7oo  lio.oo         ito.oo         iB>-:o 

IHE  Tfl 


CO. 00  90.00 

THEIR 


Ul.OO  ISO. 00 


60.00  90.00 

THEIR 


lU.DO  ISO. 00 


Figure  2.1   Differential  elastic  scattering  cross  section 
function  compared  against  experimental  data 
at  . 6eV  Seng  (1975),  and  20  and  53  eV  Trajmar 
et  al.  (1973)  . 
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Inelastic., Cross,  Sect  ions 

He   beqin  this   section  with   a   description  of   the 
differential   scatterinq   cross  sections   for   ionization 
collisions.     The  theoretical   work   on  the   ionization 
continuum  by  Massey   and  Mohr  (1933)   and   Kott  and  Massey 

(1965)   provide  a   basic  framework  for  the   calculation  of 
simple  atouic  and  laolecular   cross  sections  for  ionization- 
Green  and  Dutta  (1967)   have  developed  a  practical  method  of 
determining  the  doubly  differential   cross  sections  d  o/dS^dT 

(T  is  the   enerqy  of  the  secondary   electron)   for  complex 

molecules,  such  as  H  0,   from  sparse  experimental  data  using 

2 

the  formalism  of  generalized  oscillator  strengths  (GOS)   as 

described  by  Kott  and  Ilassey   (1965).    A  specific  analytic 

representation  of  the  GOS  for  HO   is  given  by  Green  et  al. 

2 

(1971)  and  Olivero  et  al.   (1972)   with  a  correction  to  the 
optical  oscillator  strength   part  of  the  GOS   being  made  by 
Kuccher  and   Gre^^n  (1976).     This  representation   of  the 
angular  dependence  for  ionization   scattering  is  perhaps  the 
most  complete  available  at  present;   however  because  of  its 
complexity  it  was   felt  that  a  less   ambitious  approach  for 
this  inelastic  scatterinq  type  should  bo  made.  Following  the 
work  of  ZBW  (1932)   we   will  generate  the  scattering  angles 
according  to   the  prescription   of  Grosswendt   and  waibel 

(1978).   The  functions  are  given  in  Table  2.3, 
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TABLE  2.3 
Ionization  scattering  algorithms. 


Inelastic  Scatterinq  (  >  100  eV  ) 


cos  e^  =  •l-T/E(l-E(l-T/E))/C 


cos  02  =  /1-(1-T/E)/(1-T/C)  T>200eV 

cos  92  =  .707R  T<200eV.  00362  =  -1+2R  T<SOeV 


Inelastic  Scatterinq  {  <  100  eV  ) 


cos  6   =  .707  +  .292R 


=  .707R    T>50eV 
=  -1+2R    T<50eV 
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At  low   energies  (<100  eV)    these  functions  are  rather 
arbitrary;   fortunately  at  these  enerqies  the  elastic  cross 
section  is  beginning  to  dominate   and  we  can  attribute  most 
of   the   spatial   characteristics   of   the   collisional 
distribution  to  elastic  scattering. 

While  it  is   possible  to  derive  the   differential  cross 
section  S(E,T)  from  the  GCS  a  simple  invertible  function  has 
been  developed  by  Green  and  Sawada  (1972).   It  is 


S(E.T)  =  A(E)r^(E)/([T-T  (E)]  ^+r^(E)} 

°  2.5 


where  E   and  T  are  the   incident  and   secondary  electron 
energy.    The  adiustable  parameters  r  and  Tq  have  been  taken 
from  Olivero  et  al.  (1972).  The  function  was  then  integrated 
over  T  to   give  the  total  cross  section   and  the  parameters 
Kw  K„,  J.,   and  Jq  were  adjusted  to  reproduce  experimental 
data.  Four  ionization  cross  sections  were  considered.    The 
data  for  the  1D1+3A1  cross  sections  were  taken  from  Kark  and 
Egger  (1976)    and  for  the   1B2+2A1  cross   sections  from 
Schutten  et  ai.    (1966).   The  ratios  of   the  states  were 
implied  from  the  work  of   Paretzke  and  Berger  (1978).    The 
states,    their   energy,     and   associated   analytic 
representations  are  given  in  Table  2.4.    The  functions  are 
plotted  in  Fig.  2.3  and  2.4. 
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TABLE    2.^ 
The   total  and    secondary  electron   cross   section. 


Ion.    States  Threshold  K  K 

A  B 


1B1  H20+  12.620  2.590  173.1 

3A1  H20+  14.750  2.115  176-4 

1B2  OH+,H+,0+  18.510  2.069  116.2 

2A1  H+,0+  32.400  0.884  241-0 


J 

h 

1 

.0 

1. 

.0 

1 

.0 

1 

.0 

fB I To 

IBI  H20+  -5.580        11.40        -2.37 

3A1  H20+  -10.39        11-40        -2-37 

1B2  CH+,H+,0+  -7.885        11.40        -2,37 

2A1  H+,0+  -30.00        11.40        -2.37 


Next   we   consider  the   excitation   cross   sections. 
Inelastic   scattering   cross   sections   for   excitation 
collisions  are  assumed  to  be   of  secondary  importance  to  the 
elastic  and   ionization  scattering  distributions   since  at 
hiqh  enerqy   (above  100  eV)    ionization  is   the  dominant 
collisional  process  and  below  this  energy  elastic  collisions 
doainaxe  (see  Fig.  2-2)-    The  excitation  cross  sections  are 
represented  oy   the  analytic  function 
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q  F.       W.°'i  "^3   W.  ^  j            v.. 
o.(E)=  -J-L  [l-(-=r^)    ]    C-^  )     + -^ 2.6 

where  tue  resonance   term  is  used  only   for  the  vibratioual 
states.    In  all,  fifteen  individual  states  were  considered. 
These  are  presented  in  Table   2.5.    The  functional  form  is 
taken  from  Green  and  Stolarski  (1972)   and  the  laodif icat ion 
for  the   vibrational  states  is   due  to  ZBW   (1982).    Tae 
parameterization  for  all  excitation   states  is  described  by 
Z3»  (19^2)   which  draws  on  work  by  Olivero  et  al.   (1972), 
Kutcher  and  Green   (1976),   and  data  by   Beenakker  et  al. 
(1974),   i^ohlaann  and  de  Heer   (1979),   and  Fuiita  et  al* 
(1977)  on  dissociative  excitation.    The  parameters  used  in 
Eq.   2.6  are  presented  in  Table   2.5  and  the  functions  are 
plotted  in  Fiq.   2.5.    Special  note  should  be  made  of  the 
Rydberq  states  where  contributions   from  all  possible  states 
with  n=3  or  above  are  included  according  to  the  prescription 
of  Green  and  Stolarski  (1972)  and  Jackiaaii  et  al.  (1977). 

Finally  we  calculate  the  total  inelastic  cross  section  by 
suraainq  the   aoove  functions  and  above   20  eV  we   fit  this 
cross  section  with  the  function 

^o^        E     P 
'^tOt^^>  =  ^1"  t-'^W-^^  -^3  2.7 


where  F  =  2.9b25,  :if  =  5.5343,  °'  =  . 02239,  and  P=1. 7737.    Below  20 
eV  we  enter  the  total   cross  section  discretely.    The  fit 
above   20  eV   is   virtually   indistinquishable  from   the 
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lABLS  2.5 
Paraiaeters  for  the  excited  states  of  H20. 


Exc.  States  Threshold         a  P 


(A1)  V2 

0.1980 

1.0000 

5-1930 

(Al)  VI  (B1)  V3 

0.21600 

1.0000 

6.0900 

(Bl)  V1(31)  V3 

0.4600 

1.0000 

6.0900 

TRIPLET 

4.5000 

1.0000 

3.0000 

TRIPLET 

12.5000 

1.0000 

1.0C00 

DIS.CONT 

7.4000 

1.0000 

3.0000 

DIS.CCNT 

9-6700 

1.0000 

3.0000 

DIFFUSE    BANDS 

13.3200 

1.0000 

3.0000 

H*    LYMAN    ALPHA 

17.0000 

0.6000 

6.0000 

H*    BALHEB    ALPHA 

19.0000 

0.7500 

3.0000 

OH*    3064A 

10.0000 

0.4500 

0.6500 

SESIES    A 

9.9998 

1.0000 

3-0000 

SEHIE3    B 

9.9998 

1.0000 

3.0000 

SERIES    C 

11.0600 

1.0000 

3.0000 

SERIES    D 

11.0600 

1.0000 

3.0000 

W 

Q 

F 

IaTTvI 

0.  1980 

0.2390 

0.25  4     (- 

■4) 

(Al)  V1(B1)  V3 

0.4600 

0.5000 

0.850     (- 

-3) 

(Bl) VI  (Bl) V3 

0.8990 

0.5000 

0.850     (- 

-4) 

TRIPLET 

4.5000 

3.0000 

0.7C0     (- 

-2) 

TRIPLET 

9.8100 

3.0000 

0.180     (- 

■1) 

DIS.CONT 

7.4000 

0-7500 

0.630     (- 

-1) 

DIS.CCNT 

9.6700 

0.8000 

0.  138 

DIFFUSE    BANDS 

14.7300 

0.7500 

0-626 

H*    LYi4AN    ALPHA 

17.0000 

1.3000 

0.994     (- 

-1) 

H*    BALDER    ALPHA 

18.5000 

1.0000 

0.266 

OH*    306i;A 

9.5000 

0.7250 

0.538     (- 

-1) 

SERIES    A 

9.9998 

0.7250 

0.340     (- 

-1) 

SERIES    B 

9.9980 

0.7500 

0.338     (- 

-1) 

S  zii iiXiJ  S    C 

11.0600 

0.7500 

0.771     (- 

-1) 

SERIES    D 

11.0600 

0.7500 

0.771     (- 

-1) 

Vi  V2  V3 


{A1)V2  0.135  0.677     (-1)     0.275 

(A1)V1(B1)V3  0.495     (-1)  0.130     (-1)     0-587 

(B1)V1(B1)V3  0.495     (-2)  0.130    (-1)     0.587 
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Figure  2.5   Total  excitation  cross  sections  versus 
incident  energy. 
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calculated  value.    Tte  total   elastic  and  inelastic  cross 
sections  are  graphed  in  Fiq^  2,2, 

Multiple  Elastic___Scatterinq 

Below  50  eV  the  elastic  cross   section  is  at  least  twice 
as  iarqe   as  the  inelastic  cross   section.   This   ratio 
continues  to   increase  as   the  electron   enerqy  decreases 
approachinq  the  threshold  for   the  lowest  excitational  state 
of  H2O.    The  electrons  will   then  slowly  approach  thermal 
equilibrium  with  the  H2O  molecules.  It  would  be  uneconomical 
to  follow  these  electrons  collision  by  collision  durinq  this 
staqe  of  their  degradation;    therefore  a  multiple  elastic 
scattering  distribution  (MESD)  is  used  below  20  eV, 

MESDs   have   been   used    foe   condensed   histories 
calculations,  for  example,  by  Berqer  (1963),   Jackman  (1978) 
and  Sinqhal  aud  Green  (1981).    The  latter  two  studies  used 
the  work,  by   Kutcher  and  Green  (1976),   who   carried  out  a 
Monte  Carlo  calculation  in  order   to  follow  the  position  of 
electrons   elastically   scattered   after   s   uiean   free 
pathlenqths  and  characterized  the   results  using  cumulative 
functions  for   the  radial,   lonqitudinal  and   polar  angle 
distributions. 

Following  the  work  of  Kutcher  and  Green  (1976)  (KG),   we 
have  made  a  detailed  Monte  Carlo   study  of  the  position  and 
direction  of  the  electron  distribution  after  it  has  traveled 
1  to  500  mean  free   pathlenqths,    Usinq  the  computer  code 
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developed    by      Kutciier    ve   have      essentially   solved      the    one 
velocity    Doltzmann    equation 


3f 

~(R.Z.U.A4..S)    +   Q(<t>.e)«7f  (R.z.u.A(t>.s) 

2.8 
=    [f  (R.Z.U.A4..S)    -    f  (R.z.u'  .A<t.'  .s)]    I   (In'-fil)dQ' 


where  R,  Z,  y=cos6/  A4>  =ij>  ~  <^    and  s  are  defined  in  Fiq,   A. 2 
and  oaeqa   is  a   unit  vector  alonq   the  direction   of  the 
velocity^   In  particular  we  have  calculated  the  distribution 
functions  for    E  (s)  ,   Z  (s)  ,    and  X  (s)  =  ( 1 -^  (s)  ) /2   and 
represented  them  with  analytic  functions   based  on  those  by 
KG,  with  minor  modifications,  at  nine  enerqies  from  .1  eV  to 
50  eV.     The  functional   forms  used   in  fittinq   these 
distributions  are  presented  in   Table  2.5-    Eepresentative 
qraphs  or  the  distributions  are   shown  in  Fiqs.   2.6-2. lU. 
The  solutions  qenerated   by  the  Monte  Carlo   simulation  are 
represented  as  discrete  points.    We   see  that  the  analytic 
functions   represent   these  solutions   quite   accurately, 
altiiouqh  for  fewer  than  four  nean  free  pathlecqths,  the  fits 
systematically   undershoot   or  overshoot   the   qenerated 
solutions.    This  is  of  minor  concerii,   since  these  results 
are  usually  applied  only  to  those   cases  in  which  the  number 
of  elastic   collisions  between  each  inelastic   collision  is 
much  qreater  than  tour. 
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TABLE  2.6 
Cumulative  distribution  functions  for  elastic  scatter inq. 

F(R.S)  =  1-e"*^ 

V(s)  =  2[1-Ae^^''\^  ] 

MS)  =  (-^)   .  i-^) 


F(z.s)  =  {l-e~^^[F(0)~^^''-l])-'' 


D 
-(s/s  )  ^ 
v(s)  =  1-A  e 

V 


-(s/s  ) 
F(0)  =  Fptl  -  e      "^    ] 

u(s)  =  (H  +  s^)/s'^ 


F(X.s)  =  (^  X)"X^ 


-s/s 
ir(s)  =  Le 


-s/s 
e  ( s  )  =  Me 


-(s/s  ) 
t(s)  =  1-e 


X  =  (1-cos  e)/2 


TABLE    2.7 
Parameter    values    for    radial   distribution. 
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E(eV) 

.2 

.5 


B 


5.24 
5.37 

4.74 


Y 

1.00 
1.  10 


.740 
>728 
.739 


2.79 
2.87 
2.73 


10.0         7.09  1.G9 

20.0         6.88  1.11 

50.0         9.00         1.07 


W 

TToo" 

1.95 

2.06 


1.0  5.94  1.00  .641  2.94         2-21 

2.0  7.00  1.00         .596         3.09         2.28 

5.0  6.37  1-02  .664         2.36         2.07 


.623         3-05         2.07 
.645         3.09        2.00 


TToo" 

.  100 
.100 


w 


.734 
.755 
.755 


-100  .790 

.100  .786 

.100         ,760 


.63 


3.73         2.32 


-100 
.100 

.  100 


-753 
.731 

.7  99 
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TABLE    2.3 
Parameter    values   for   loiiqitudinal    distribution- 


E(eV) 

.2 
-5 


D 


H 


V  V  V 

'1 0  75  "173  4  ^7  50 ITT  5"50  795  7"~T  ~  5  o" 

9.90  1.50  -.252    .600  1.00       1.80 

9.20  1.40  .0085    .650  1.05      2-30 


^F 


B. 


F  "F         "F 

r435~~9722"7600' 

.435       11.5    . 640 
.470       13.0    .720 


1.0  9.15       1.52       .170       .680       1.12      2.40       .447       16.5    .719 

2.0  8.40       1.71       .350       .700       1.17       2.30       .435       16.1    .725 

5.0  8.00       1-87       .400       .740       1-17       2-60       -440       12-3    -704 


10.0  7-50       1.94       .451       .800       1.22      2.50       .428       10.2    .700 

20.0  8.00       1.75       .500       .900       1.29       2.60       .435       9.10    .750 

50.0  9.00       1.72       .475       .940       1.39       2.60       -425       14.4    .750 
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TABLE  2.9 
Parameter  values  for  Polar  Anqles  distribution. 


E(eV) 


1 


.2 

.5 


'73  oo" 
.700 
.580 


L 

""550" 

,800 
1.  10 


c 

'TasT 

>600 
.840 


M 

"IToo" 

1.60 
1.40 


1 .  oT 

1.46 
1.76 


N 

T6  50' 

1.07 
1-25 


1.0 
2.0 

5.0 


.  523 
.459 

.500 


1.28 
1.26 
1.25 


1.13  1-13  2.58  1.27 
1.50  1.10  2.84  1.08 
1.30     1.30     1.98     1.23 


10.0 
20.0 
5  0.  0 


.540 
.770 
1.40 


1.20 
.973 

.550 


1.10 
.  315 

.550 


1.40 
1.68 
2.00 


1.85 
1.60 


1.20 


1.20 


1.  19 
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P4 


X  =  (l-COS9)/2 


Figure  2.12  Cumulative  polar  angle  distributions  for 
multiply  scattered  electrons  with  energy 
.2eV  for  various  pathlengths,  S. 
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Figure  2.13  Cumulative  polar  angle  distributions  for 
multiply  scattered  electrons  with  energy 
2 . OeV  for  various  pathlengths,  S. 


1.0 


43 


1.0 


.8  - 


6  - 


Cm 


4  - 


2  _ 


.4         .6 

X  -  (i-cose)/2 


1.0 


Figure  2.14 


Cumulative  polar  angle  distributions  for 
multiply  scattered  electrons  with  energy 
20eV  for  various  pathlengths,  S. 
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Xkj§_lisi:^_S£ectrum_aiid_Yieids 
Tiie   concept   of    a    yield      spectrum    was   introduced    by    Green, 
Garvey    and   Jackinan    (1977)        (GGJ)  .       This   concept   traces    its 
evolution   from    the    work   on      the   equilibrium   £lux    of    Spencer 
and   Fano     (1954)^    Peterson    (1969)     and    Jura    (197  1).       The   yield 
spectrum   has      properties    whicii   are      much   simpler      than   tie 
equilibrium  electron      distribution.        Its      shape   can      be 
described    by   a    continuously    varyinq    function    of    incident    and 
spectral    electron    enerqy   and   as      demonstrated    by   GGJ     (1977) 
this   functional   dependence      is    similar   for    many      atoias   and 
molecules.         As   introduced    by   GGJ    (1977)     the    yield    spectrum 
is   defined   by 

U(E.Eq)     =    N(E.E^)/AE    =    U^(E.E^)H(E^-E-Ejj^)+«(E^-E)  2.9 


where  N  is   the  total  number  of   inelastic  collisions  which 
occur  in  the  enerqy  interval  A3  centered   at  E  and  E   is  the 
enerqy  of   the  incident  electron,     II  is   the  Heavisida 
function  with  E^       set  at  the  lowest  threshold   value  of  all 
the  inelastic  interactions   considered  and  *  is   the  Dirac 
delta  function  which  allows  the  second  term  to  represent  the 
source  contribution  to  the  yield  spectrua. 

This   definition  has   also  been   qeneralizcd  to   four 
dimensions  (GS) .   The  functional  form  is 
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U(E.p.X.E  )  =  N(E.p,z)/AE&sAz 

=  U^(E.p.z.Eq)H(E^-E-Ejjj)+D(p.z.E^)«{E^-E) 


2.10 


where 

and  N(E,  p,  z)   is  the  total   nuniher  of  inelastic  collisions 

which  occur  iii  the  volume  element   As,  z  centered  at  p,  z  and 

the  enerqy  interval  Ae  centered  at  E-    Analytic  expressions 

for  U   and  D  are 
a 

U^(E.p.Z.E^)  =  I    A.(E^)G.(p.z.E^Xi(E.E  JD^^ 


DCp.Z.E^)  =  A4(E^)G4(p.Z.E^)D^^ 


2.11 


waere 


A.  (E^)  =  {D^  j"  r   G.  (p,z.E  )pdpdz) 


Xi(E.E^)  =  Ci-^  X.E,e 


\^^i 


and  D   is  the  density  of  the   mediua  with  units  of  qrams  per 
o 

centimeter  cubed.    Note,   by   dcfininq  the  normalization 

constant   A  with   D"^    explicitly   displayed,   we   have 

o 

essentially  taken  all  density  dependence  out  of  it  since  the 
inteqrated  quantity  is  proportional   to  the  inverse  density 
cubed.   The  paraaeterization  for  li^O  is  qiven  in  Tables  2.10 
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and  2.11-     Representative  qcaphs  of   the  two   and  four 
dimensional  yield  spectra  are  given  in  Fiqs.  2.  15-2.  17. 

It  is  also  necessary  to  consider  the  ultimate  fate  of  the 
electrons  both  the  incident  and  those  produced  in  ionization 
collisions  in   order  to   carry  out  the   next  part   of  the 
calculation.    The  MC   code  was  modified  to   follow  these 
electrons  below   the  lowest   excitation  threshold   enerqy. 
These  electrons  undergo  elastic  collisions,  but  it  is  a  well 
known  result  froia  classical  mechanics  that  they  can  exchange 
kinetic  enerqy  with  the  B  JJ   molecule,   losing  enerqy  on  the 
average,  until  they  reach  thermal  equilibrium.    Because  the 
average  enerqy  lost   to  the  H  C  molecule   by  tne  electrons 
approaches  zero  as  the  electron  approaches  equilibrium,   it 
is  not  practical   to  follow  the  electron   to  equilibrium  iu 
the  MC  simulation.    Therefore  a  cutoff   or  E  =  .125  eV  was 
chosen;   that  is,  once  the  electrons  enerqy  fell  below  this 
energy  its   position  was   recorded  and   it  was   no  longer 
followed.    A  new  guantity  is  defined,  called  tiie  deposition 
spectrum  which  is   defined  as  the  number   of  electrons  with 
energy  less  than  .125  which  are   in  the  volume  element  As,  z 
centered  at  p,  z.    This  deposition   spectrum  is  unlike  the 
yield  spectrum  which  represents  the   number  of  collisions  or 
the  deposition  spectrum  as  defined  by  Singhal  et  al.   (1983) 
which  represents  the  number  of  electrons  with  energy  E  after 
their  last  inelastic  collision.   The  deposition  spectrum  can 
be  analo^ou3ly  represented   by  a  function  similar   to  that 
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lABLE  2. 10 

Parameters  for  two  dimensional  yield  spectruia. 


U(E.E    )     =    I     x:   (E.E^) 
i  =  l  ° 


Xi(E.E^)    =   Ci.   \^e"''°'^i 


E.     =    E^/1000 

Iv  O 


^i  ^i  ? 


1 


1  -0339  .8036  40.3 

2  0.0  99-85  4,61 

3  0.0  830.6  .074 
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TABLE    2. 11 
Parameters    for   four    dimensional    yield    spectrum. 


G^(p.z.E^)    =    e 


-B  .D  /p  \(>   .(z-z   .   )   ^ 
'  1    o  ^  1         oi 


p.(E^)     =    A./{E^-.B.)ln(E^-HC.) 
z^.(E)    =    (E^.B.)/A.ln-\E^.C.) 


1 
2 
3 

'T 

2 
3 

4 


A 


21-i;3  (7) 
20.  32(7) 
13.30(7) 
91.92(7) 
Il573l7r 
619.6  (7) 
34.07(7) 
88.70  (7) 


B 


-75.90 
-25.80 
-5-224 
87.39 

"looTs 

465.1 

-52.95 

355-6 


C 

-92.63 
-90-97 

6.  005 
'l0574 
833-  1 
8.82 
-66.04 


1 
2 
3 

4 


T7091' 
.5902 
.  2835 
.700 
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Figure  2.15 


E(eV) 

Two  dimensional  yield  spectrum  and  electron 
deposition  spectrum  as  function  of  spectral 
energy  E  at  various  incident  energies,  £„ . 
The  points  represent  the  distribution  values 
obtained  directly  from  the  Monte  Carlo 
Transport  code. 
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Figure  2.16   Four  dimensional  yield  spectrum  for 
electron  with  incident  energy  500eV 
versus  spectral  energy  E  at  the  spatial 
points  (p,  Z)  . 
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Figure  2.17 


Four  dimensional  yield  spectrum  for  electron 
with  incident  energy  200eV  versus  spectral 
energy  E  at  the  spatial  points  (p,  Z) . 
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used  for  the  four  dimensional  yield  spectrum  except  that  it 
requires  only  one  term  instead  of  three  and  has  no  spectral 
energy  dependence.   The  function  is 


DS(p.z.E^)  =  N(E^)A(E^)G(p.Z.E^) 


where  the  incident  enerqy  dependence  is  qiven  by 


2.12 


GCp.z.Eq)  =  e  2.13 


and  the  parameters  are  qiven  in  Table  2.12.    Normalization 
factors  A   for  yield  and   deposition  spectra  are   qiveu  in 
Table   2.13.    Tue   function   is   plotted  for   various 
isoelectronic  densities  in  Piq.  2.18  and  Fiq.  2.19. 

Havinq  obtained  the  yield  spectrum   we  can  then  calculate 
the  specific  state   yields  {#iflol.   of  H  0   in  a  particular 
state  per  cubic   centimeter  due  to  all   collisions  for  one 
electron)  as  a  function  of  position  by  inteqratinq  the  yield 
spectrum   multiplied  by  tae    probability   of  a   collision 
producinq  a   molecule  in   the  i^^  state   over  E   troai  the 
threshold  of  that  particular  state   to  the  incident  electron 


enerqy.   Analytically,  the  yields  are 

E 
°  P^(E)U(E.p.z.^  )dE  2.14 

i 


J.(p.z.E^)  =  I^°  p.(E)U(E.p.z.^  ). 
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TABLE  2. 12 
Parameters  for  electron  deposition  spectrum. 


-PD  /p^+i{Z-2^)  ^ 

G(p.z.E^)  =  e    °         ° 
P(E^)  =  A/(E^+B)ln(E  +C) 

-(E  /E  )" 
6(E^)  =  1-Ae    °   ^ 

^o^^o^  =  (E^+B)/A  ln-^(E^+C) 


A=19.30(7)    B=-22.77       C=-91.70 


^o  A=102,8(7)    B=213.7        C=-90.53 


A-,3809(7)    E=5n.9        N=a,98 


Ene  rqy 
(eV) 

100 

200 

500 

1000 

Number  of  thermal  electrons 
produced  per  incident  electron 

3.81 

7.00 

16.62 

32.87 
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TABLE    2.  13 
Noriaalization   constants,    A- 


Incident  lield    Spectrum  Electron 

Enerqy     (eV) 

12  3  4  Spectrum 


100  3.25    (19)     -399    (19)     .100    (19)     3.35     (20)     .3774(19) 

200  1.49     (19)     ,412     (19)     .110    (19)      1.39    (20)     .3400(19) 

500  ,0796(19)     .0434(19)     .0140(19)     -2590(20)     .4927(13) 


1000  .112    (18)     -729     (17)     .250    (17)     .5587(19)     -9711(17) 

2000  .169     (17)     .120     (17)     .422    (16)      .1051(19)      .1597(17) 

5000  .144     (16)     -106     (16)     -382     (15)     .1020(18)     -1420(15) 
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where  J  •  is  the  yield  for  the  i    state,  H  ^  is  the  threshold 
for  the  i   state  and  p^  is  tne  probability  of  interaction 
for  this  state,  qiven  by 


p^(E)  =  <^i(E)/o^Q^(E)  2.15 

Pert oriainc,'  the  integration  over  all   space  as  well  qives  the 
number  of  collisions   of  a  particular  kind   which  occur  for 
each  electron  as  it  is  degraded.    Of  course  in   this  case  it 
is  simply  easier   to  start  with  the   two  dimensional  yield 
spectrum.    Reqardless,  it  is  then  possible  to  calculate  W, 
the  energy  loss  per  ion  pair  produced.   We  have  calculated  a 
value   of  W(9V/ip)    for  H  O   which   is  comparable   to 
experimental  rasuits.    The  results  are  presented  in  caapter 
IV.   A  few  representative  isoyield  contours  are  presented  in 
Fiqs.   2.20  -2.23.    Ihese  will  te   used  as  input  into  the 
diffusion  kinetic   calculations  which  follow.     In  these 
calculations  more   than  one  state   has  contributed   to  the 
production  of  a   specific  species.    The  states   and  their 
products  are  presented  in  Table  2.14.    It  is  also  possible 
to  calculate   isodose  curves  usinq   the  yielu   spectrum  by 
perioriuinq  the  followinq  sum 


Dose{p.z.E^)a  I   W.  J.  (p.z.E^) 


2.16 


over  all  states. 
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TABLE  2-14 
States  and  associated  products. 


states 
7aT7v2 

{^^)  VI (Bi) V3 

(Bl) VI  (BI) V3 

TRIPLET 
TRIPLET 
DIS^CONT 
DIS.CONT 
DIFFUSE  BANDS 
H*  LYHAN  ALPHA 
H*  BALtlER  ALPHA 
OK*  3064A 
SERIES  A 
SERIES  B 
SERIES  C 
SERIES  D 


products 

> 

M20* 

H20* 

H2C* 

H20* 

H20* 

H,     OH 

H2.    0 

H20+, 

e- 

H,  on 

H,    OH 

H,    OH 

H20  +  , 

e- 

H20+, 

e- 

a2o+. 

e- 

fi20+. 

e- 

1B1  B2  0+,  e- 

3A1  li2G+,  e- 

132  OH+,H/  H+,CH/  0+,H2 

2A1  H+,OH/  0+,H2 
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CHAPTER  III 
DIFFDSION  AND  KINETICS 


General 
After  calculatinq  the  yields  of   tiie  various  excited  and 
ionized  states  of   H2O/   which  occur  in   less  than  10 (-13) 
seconds,   we   are  interested  in   studyinq  the   effects  of 
diffusion  and  kinetics  on  these  chemicals.    Usinq  standard 
diffusion  theory   it  is  possible   to  calculate   the  steady 
state  distribution  of  the  initial  species  and  those  products 
like  H2O2   which  are  created   by  the  interaction   of  these 
initial  products.    In   order  to  accomplish  this   we  have 
developed  a  proqraai   called  DIFKT  which  can   calculate  time 
dependent  and  steady  state  spatial  concentration  profiles  of 
all  the  initial  species  and  their  products.    Details  of  the 
proqraui  are  addressed  in  chapter  III. 

The  use  of  a  diffusion  kinetic  (DK)   model  to  explain  the 
dependence   of  final   products   on   the  initial   spatial 
distribution  of  radiation  products   was  first  formulated  by 
Jaffe  (1913,1929)  and  Lea  (194b).    Samuel  and  Kaqee  (1953) 
extended   this   work,    applyinq   several   mathematical 
approximations  to  laake  the   problea  analytically  tractable. 
Waqee  (1955)   considered  the  effects   of  dose  rate  for  this 
aodfel-    That  is  the  case  when   the  enerqy  is  absorbed  at  a 
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en 

sufficiently  aiyh  rate  that  we   cannot  treat  the  problen  as 
the  sum  of  the  effects  of  individual  particles.    Kuppermann 
(196  1)   treating  the   problem  in  the  absence   of  dose  rate 
effects  removed  the  mathematical  restrictions  imposed  in  the 
previous  work.    He   solved  a  system  of   diffusion  kinetic 
equations  of  the  form 


ac.(r.t)=D.72c.(r.t)-Z  C^(r . t)C^(r . t)k^^.  3.1 
nm 


3t 


C^Cr.O)  =  C^^Cr) 


where    D.,         k.  and  C-       are    respectively      the   diffusion 

i'  nmi  1 

th  -  ^       .  th 

constant    of    the    i         species,      tne    rats    constants    tor    the    i 

th 
species    and      the    concentration      or    the      x         species      as    a 

function    of    time      and   space    (Kuppermann,       1959,  1961    ; 

Kuppermann   and   Bolford,     1962a,     1962b) ,         The   initial    species 

present   and      their    spatial      concentration    were      considered 

somewhat    arbitrary.  Recently    this      approach    has      been 

extt^nded    by    Ilaqee    and   Chatterjee     (1978,1980)       and    Chatteriee 

and    [laqee    (1980).       They   essentially   start    their    calculations 

after    ail   ions      and    excited    states   have      reacted    and    formed 

neutral    qrouud    state    molecules.       The    approximations   involved 

in    these    calculations   are      numerous;       the    most   questionable 

are    those    iavolvinq    the   spatial      distribution    of    the    primary 

neutral    species.       They    are    adiusted    to   qive    the    best    results 

in      reproducing   experimental      data      for    differential      and 
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integral  G{ip/eV)   values.    There   are  no  guarantees  that 
their  choices   for  the  primary  distributions   represent  the 
correct  input  for  the  problem,   though  the  results  presented 
for  a  Fricke  dosimeter  system  compare  well  with  experiments. 
In  general  this   approach  does  not  address   the  early  time 
chemistry. 

A  more  basic  approach  to   the  problem  of  calculating  the 
initial   species  present   using  Monte   Carlo  methods   was 
initiated  in  the  late  70 's  at   Oak  Ridge  National  Labs  (Hamm 
et  al.,  1978;  Ritchie  et  al. ,  1978)  and   Los  Alamos  National 
Labs  (ZDW,  1982;   Zaider  and  Brenner,   1983  and  Brenner  and 
Zaider,  1983).    This  approach  was  recently  connected  to  the 
work  of  Magee  and  Caatterjee   (1980)   (Wrignt  et  al.  ,   in 
press;   Wriqht   et  al.    1983).    Initial   prediction  of 
experimental  G  values  are  reasonable. 

While  a  tine  dependent  solution  for  one  electron  track  as 
presented  is  very  useful  in  predicting  experimental  results, 
a  continuous  beam   (or  at  least  pulsed)    of  electrons  is 
normally  encountered.     The  raaior   part  of   the  present 
analysis  has  concentrated  on  looking  at  a  continuous  beam  of 
electrons.    The  transient  behavior  of   this  system  has  not 
been  addressed  in  the  diffusion  kinetic  calculation. 
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Kinetics 

In  solving  a   system  which  involves  both   diffusion  and 
kinetics  it  is  often   instructive  to   solve  the   problem 
iqnorinq  diffusion.   In  this  section  we  present  the  chemical 
system  which  exists  iiaiaediately  after  the  incident  electrons 
have   deposited   their  energy.     This   system   exists 
approximately  10  (-13)   seconds  after  an  electron  is  iuiected 
into  the  medium.    Thus  we  see  that  the  physical  part  of  the 
deposition  process   finishes  very   quickly,    A   complete 
description  of  these  events  is  presented  iii  chapter  I. 

The   system   of  nonlinear   first   order   differential 
equations  which  describes  the  chemistry  is 

^^i^^^  =  I   k   .C  (t)C  (t)  ,  7 

i; ^    nmi  n^    m^  '  3.2 

dt        nm 

C.  (O)  =  C  . 
1  ^  "^      01 

Table  3^1  lists   the  species  considered  and   Tables  3.2  and 
3-3   list   the   reactions.    The   system   consists   of 
electrons, ions  and  excited   species  of  1-12  0  as   well  as  the 
dissociated  products.    The  neutral  chemistry  reaction  rates 
have  been   taken  from  a   recent  coicpilation   by  Westbrook 
(1931) .    The  ion   neutral  reaction  rates  are   taken  from 
Huntress  (1973)   and  Ferguson  (1973).    The  ion  clustering 
reactions  are  taken  from  the  work   of  Good  et  al.   (1970). 
The  ion  recoabination  rates  froia  Leu  et  al.   (1973).    The 
chemistry   is  fairly   well  known   for  H2 0   vapor  and   is 
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delineated  in  the  books  by   Allen  (1961),   Spinks  and  Woods 
(1976),  Venuqopolan  and  Jones  (1968)   and  articles  by  anbar 
(1968),  Thomas  (1969)   and  Hunt  (1976).    The  most  serious 
deficiency  is  the   lack  of  information  on  the   fate  of  the 
excited  states  of  112  0.    In  a  recent  study  for  liquid  water 
these  states  were   assumed  to  dissociate  in   the  same  tioe 
scale  as  the   interaction  of  the  H2O+  ion   (Wriqht  et  al . , 
1933)  .    For   vapor  we   would  expect   a  sliqhtly   slower 
dissociation  rate.    However  the  information  ia  this  area  is 
hiqhly  deficient  (Inokuti,  1983). 

The  system  of  differential  equations   is  solved  usinq  the 
Gear  method  (Hindmarsh,  1974).    The  proqram  is  based  on  the 
routine  GEAR  (IM3L,  1981)   and  includes  an  adaptation  of  the 
kinetic  subroutines  from  the  diffusion  kinetic  proqram  DIFKT 
(see  chapter  III)   in  order  to  set  ap  the  kinetic  equations. 
Usiiiq  the  results  froia  the  calculation   of  yields  to  set  the 
relative  concentrations   (see  Table  3.1)   of   the  species 
initially  produced  by  the  incident  electron,   it  is  possible 
to  study  the  time  development  of  these  concentrations.    The 
relative  concentrations   are  based  on  the   distribution  of 
products  presented  in   Table  2- 14  and  the   thermal  electron 
distribution.   The  reaction  mechanism  is  presented  in  Tables 
3.2  and  3.3.    The  results  of  tae  calculation  are  presented 
in  Fiq.  3.1  and  discussed  in  chapter  IV. 
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TABLE  3. 1 
Chemical  species  included  in  the  diffusion  kinetic  system. 


Species 


Diffusion 

Coefficients 

{cm**2/sec) 


E 

5.05 

H2 

.8644 

02 

.2305 

H202 

.2327 

H02 

.2339 

OH 

.3695 

0 

.3754 

H 

1.292 

H+(n20)3 

.070 

U20  + 

.250 



H  + 

1.290 

H+  (H20)  2 

.  1  10 

H30  + 

.245 

0H  + 

,370 

;i20 

.250 

Initial  Concentrations 
Kinetic  Systeni 
(moles/cm* *3) 


9.900 

2.643 

0.0 

0.0 

0.0 


1.395 
6.044 
6.215 
1.572 
1.101 


1.091 
1.581 

0.0 
0.0 
0.0 


(-08) 
(-12) 


(-12) 
(-12) 
(-12) 
(-11) 
(-11) 


(-12) 
(-12) 
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lABLE    3.  2 
Reaction   rates   at   temperture   equal   to   400    K. 


Reaction  Reaction    Equation 

Number 

1  H    +    02    ->    0    +    OH 

2  0    +   OH    ->    H    +    02 

3  H2    +    0    ->    H    +    OH 

4  H    +   OH    ->    H2    +    0 

5  H20    +    0    ->    OH    +    OH 

6  OH    +    OH    ->    H2C    +    0 

7  H20    +    H    ->    H2    +    OH 

8  H2    +    OK    ->    H20    +    H 

9  H2C2    +    OH    ->   H20    +    H02 

10  H20    +    H02    ->    H202    +    OH 

1 1  H20    +     H    ->    H    +    OH    +    M 

12  H    +    OH    +    ;i    ->    H20    +    a 

13  H    +    02    +    M    ->    H02    +    M 

14  H02    +    H    ->    H    +    02    +    M 

15  H02    +    0    ->    OH    +    02 

16  OH    +    02    ->    H02    +    0 

17  H02    +    H    ->    OH    +    OH 

18  OH    +    OH    ->    H02     +    H 

19  H02    +    H    ->    H2    +    02 

20  H2    +    02    ->    H02    +    H 

21  H02    +    OH    ->    H20    +    02 

22  H20    +   02    ->    H02    +    OH 


Reaction    Sate 
cm, moles, sec 


1.244  (    5) 

.6294  (13) 

.9954  (    8) 

.5300  (    9) 

3.187  (    3) 

.7910  (12) 

7.729  (    2) 

.3359  (11) 

.  1038  (13) 

.3432  (-4) 

.9265  (-41) 

.8812  (18) 

.5837  (16) 

,  1902  (-9) 

.1424  (14) 

.7544  (-17) 

.2299  (14) 

.1472  (-08) 

.  1041  (14) 

.1434  (-17) 

.  1424  (14) 

.2760  (-25) 
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TABLE    3.3 
Reaction    rates   at    temperture   equal   to   400° K. 


Reaction  reaction  Equation  UToleVslT 

lluDiber 


23 
24 
25 
26 
27 


H2C2    +    02   ->    H02    +    H02  .2000  (-9) 

H02    +    H02    ->    H202    +    02  ,2843  (13) 

H202    +    tl    ->    OH    +   on    +    M  .1651  (-7) 

OH    +    OH    +    M    ->    H202    +    M  ,5380  (18) 

H2C2    +    li   ->    H02    +    112  .1514  (11) 

28  HC2    +    H2    ->    H202    +    H  43.84 

29  0    +    H    +    E   ->    OH    +    M  I'OOO  (16) 

30  OU    +    il    ->    0    +    H    +    M  .4191  (-39) 

31  02    <■    a    ->    0    +    0    +    H  .7448  (-47) 
0    +   0    +    K    ->   02    +    :i  .8743  (15) 


32 

33  H2    +    11    ->    fi    +    H    +    M 


7659       (-38) 


34  H    +    K    +    M    ->    H2    +    a  .30  20       (16) 

35  1120+    +    n20    ->    H30+    +    OH  .1130       (16) 

36  0H+    +    H20    ->    H20+    +    OK  .3480       (15) 

37  0H+    +    n20    ->    030+    +    0 

38  n+    +    n20    ->    H20+    +    H 

39  H30+    +    H2C    +    M    ->      +H(n20)2    +    M 

40  +H(H20)2    +    K20     +    K    ->       +H(H20)3    +    « 
i;l  +H(H20)2    +    S >    H    +    H20    +    H20 

42  +H(H20)3    +    K >    H    +    020    +    H20    +    H20    .2770 

43  H20+    +    H2    ->    H30+    +    H  .3670       (15) 
1^4  H30+    +    2-    ->   H20    +    H  .7830       (18) 


.2839 

(15) 

.4939 

(16) 

.  1230 

(22) 

.8339 

(21) 

.  1630 

(19) 

.2770 

(19) 
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i, I  J  Mini:   4  1,  I  i I  I  ir 4  iiMiiii  4111 )^ !,"n-^.""'?o-^-""'i"n-'-"'"Trf  • 


REfiCTIQN  TIME,  SEC 


Figure  3.1   Results  from  kinetics  program  concentrations 
versus  time.   Based  on  system  presented  in 
Tables  3.2  and  3.3. 
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Develq2aent_of_the_Dif fusiou_KiuGtic_SYS 
The  DK  system  as  formulated  in  this  study  is  described  by 
a  system  of   n  nonlinear  second  order   partial  differentia]. 
equations  (Green  and  Rio,   1983)  .    The  yields  discussed  in 
the  previous   chapter  are   converted  to   yield  rates   by 
multiplying  them  by  a  rate  term  (*elec./s).    As  previously 
stated,    the  incident   electron  qives   rise  to   spatial 
distributions   of  H20+,CH+,0+   and   H+   alonq  with   the 
associated    dissociated   products    and    electrons. 
Additionally,  various  excited  states  of  H2O  are  produced. 

The  theraal   electrons  consist  of  the   source  electrons 
plus  those   electrons  which  are   produced  in   the  ionizing 
process  and   are  degraded   below  .125   eV.    The   source 
electrons,  of  course,   have  no  companion  ion  as  exhibited  by 
the  followiiiq  equation: 

00  00  ]£ 

J„  J       [J      '^DS(p.z.E.E^)dE    + 

00  3.3 

^o 

1  J         p.  (E)U(E./D.z.E    )dE]pdpdx.^l 

i  =  l  1  ° 


The    first      term   is      the   number      of    thermal      electrons    per 
incident    electron,         the    second    terra    represents      the    total 
number    of    ions    created    and    the      third    term    accounts    foe    the 
source    electron.       Since    the    system    of    differential    equation; 
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includes  an   equatiou    for   the    electrons,      we   need   to   consider 
what   the    ultimate      fate  of   these    source      electrons   will   be. 
Dltimately,       they      will   encounter    a      qround    such      as   the 


c 


ontaiument  vessel   for  the  ^'2^'         ^^      ^^    possible   to 


hypothesize  a   reaction  which  eliminates   these  electrons, 
such  as   electron  attachment   to  H2O   (Hanrahan,   1983). 
Alternately  we   may  simply  iqnore   these  electrons  in  the 
deposition  spectrum.    For  a  1000   eV  electron  this  is  the 
same  as  reiaovinq  one  electron  out  of  thirty.    At  most  this 
would  slow  up  the  ion  chemistry  sliqhtly. 

The  physical  system  is  described   by  the  followinq  set  of 
equations: 


3C 


^(p.z.t)  -  D^y  ^  ^(p.z.t) 


3t  3,4 

-^   ^nmi^n^^'^-^^^m^^-^'^^  "^  Y.(p.z.t) 
n ,  m 

where  the  initial  conditions  are 


C^(p.z.O)  =  C^^(p.z) 


Yj^(p.z.t)  =  J^(p.z)  Rate  (tte  /sec) 

In  the  steady  state  case,  'X:/   %  is  set  equal  to  zero.    The 

required  boundary   conditions  in  this   case  are  that  the 

concentrations  equal  zero  at  infinity  and  their  derivatives, 

with  respect  to  P,       at  P    equal  zero  be  equal  to  zero.    The 
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symbol  Y.  is   used  to  denote  tue  yield  of   the  i'   species. 
The  remaininq  symbols  have   been  previously  defined.    The 
species  involved  are   listed  in  Table  3.1,    the  reactions 
considered   are  listed   in   Table   3.2,   the   diffusion 
coefficients  are  presented  in   Table  3.  1  and  representative 
yields  are  presented  in   Fiqs.   2.18-2.23.    The  diffusion 
coefficients  for   the  molecules   and  atoms   are  based   on 
qaseous  diffusion  theory  as  presented  by  Hirschfeider  et  al. 
(1954)  and  Monchick  and  Mason  (1961)   and  experimental  data 
compiled  by  Marrero  and  Mason  (1972),    The  coefficients  are 
calculated  according  to  Chapinan-Enslcoq  tiieory ,   implementinq 
the  molecular  paraiaeters   given  in  the  paper  by   Kee  et  al. 
(1933).   The  diffusion  coefficient  of  the  electron  in  {l20  is 
based  on  the   total  elastic  cross  section   as  presented  in 
chapter  II.     The  complete  set   of  parameters   needed  to 
represent  the  yields  in  the   diffusion  "kinetics  program  are 
presented  in  Tables  2.11,  2.12  and  2.13. 

In  order  to  make  the  problem  laore  tractable  the  variables 
have  been  transformed  to   dimensionless  coordinates.    The 
transformations  are 


P'=p/ap^  .z'=z/bz^  .  D.'=D./D^  ^and  C.-=C./bC         35 

2  ^ 


The  transiorined  equations  and   their  associated  initial  and 
boundary  conditions  are 
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3C'  (p'  .z'  .t)=DW'  ^C'^  (p'  .2'  .t)  +YJ(p'.Z') 


3t 


3.6 


■^^   '^;nni  C;^(p'.z'.t)C;^(p'.z'.t)  =  0 
nm 


,■'.  .',y 


2       2. 
a  p   b 

k-    -     °    k   . 
nmi    D-,  -,   nmi 

where  Y   is  qiven  by 
i 


2   2 
a  p 

YMp'.z')  =   —5 Y.(ap  p.bz  z)  3.7 

-H,0-H,0 

In      the   steady      state   case      we      make    the      additional 
t  rail  sf  o  riaa  tiofl 


2         -1 
x=  -tan     p'  3.8 

ir 

2^       -1     , 
Y=   —tan      z' 

■•  IT 

which   reduces  the   domain   of   the  problem   to   finite 
dimensions.    The  equations  can  be  restated  in  the  followinq 

form: 

3C  ' 

^     Ci    ~    I         k'    .CL  C'_      +  Yi  =  0  3.9 

n.m 


^  ^  \x.y)^i  -  ^,  „^'nmi^n   "m     "i 


:.'   =  Ci  (x.y.t)  Y|  =  Y^Cx.y.t) 


76 


wnere 


and 


„2  4  2,  irx.  a        ,  2  ,irx  .3      , 

^x.y)  =  "1^°^  ^i~>-i5^  f^°^  ^-2^-^^  * 


^^!^,^     ,^--^;,^^     1^     .  ^COS^  (^  )  |-   [COS  2(  f  A] 
irsin    {irX/2)    dX  ir  Z       ay  ^    ay 


2      2,,  2      2 
S    =    a    p^   /b    z^ 
o  o 


The  steady  iitate  solution  is   valid  for  restricted  values 
of  the  yields.   This  is  discussed  in  cliapter  IV.    Solutions 
for  the   steady  state   equations  are   presented  in   Fiqs. 
3.2-3.14,    The  sif,iulation   was  made  at  a   temperature  of 
400  °K-    The  density   of  the  medium  was   set  at  1.795  (-7) 
molecules  per  centimeter  cubed  and  a,  h    and  p   were  set  to 
1,  1,   and  1.45  respectively.    The  computational  time  was 
approximately   two  hours   and   involved  ten   iterations. 
Converqence  for  ten  iterations  was  nominal,  however  the  cost 
involved  in   extending  these  calculations   was  prohibitive. 
In  the  future,    modifications  to  the  DIFKT   proq'ram  will 
enable  it  to  make  use  of   the  knowledqe  acquired  in  previous 
calculations  and  reduce  the  computational  cost. 
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Figure  3.2 


Concentration  profiles  {p   =  0)  of  the 
electron  distribution  obtained  from  the 
steady  state  solution  of  the  diffusion 
kinetics  system.   The  solutions  are 
normalized  by  Eq .  3.5. 
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Figure  3.3 


Concentration  and  yield  rate  profiles  (z  =  0) 
(dashed  line)  of  the  electron.   Yields 
normalized  to  maximum  value  and 
concentrations  normalized  to  (0,  0)  point 
value . 
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Figure  3.4   Concentration  profiles  (p  =  0)  of  the 

H  O"*"  distribution  obtained  from  the  steady- 
state  solution  of  the  diffusion  kinetics 
system.   The  solutions  are  normalized  by  Eq 
3.5. 
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Figure  3.5   Concentration  and  yield  rate  profiles  (z  ^  0) 
(dashed  line)  of  HO"*".   Yields  normalized  to 
maximum  value  and  concentrations  normalized 
to  (0.  0)  point  value. 
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Figure  3.6 


Concentration  profiles  (p  =  0)  of  the 
H  O"^.  H'^CH  0)„  and  H'*"(H  O)   distributions 
oBtained  fcom  the  steady  state  solution  of 
the  diffusion  kinetics  system.   The  solutions 
are  normalized  by  Eq.  3.5. 
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Figure  3.7 


Concentration  profiles  (p  =  0)  of  the 
atomic  hydrogen  distribution  obtained  from 
the  steady  state  solution  of  the  diffusion 
kinetics  system.   The  solutions  are 
normalized  by  Eq .  3.5. 
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Figure  3.8 


Concentration  and  yield  rate  profiles  (z  = 
of  atomic  hydrogen.   Yields  normalized  to 
maximum  value  and  concentrations  normalized 
to  (0,  0)  point  value. 
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Figure    3.9 


Concentration  profiles  (p  =  0)  of  the  OH 
distribution  obtained  from  the  steady  state 
solution  of  the  diffusion  kinetics  system. 
The  solutions  ace  normalized  by  Eq.  3.5. 
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figure  3.10 


Concentration  and  yield  rate  profiles  (z  =  0) 
(dashed  line)  of  OH.   Yields  normalized  to 
maximum  value  and  concentrations  normalized 
to  (0,  0)  point  value. 
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Figure  3.11 


Concentcation  profiles  (p  =  0)  of  the 
atomic  oxygen  distribution  obtained  from  the 
steady  state  solution  of  the  diffusion 
kinetics  system.   The  solutions  are 
normalized  by  Eq .  3.5. 


87 


10 


-2 


c: 

•O 
(I) 
N 

•H 

i-< 
(T3 

e 

U 

o 

c 


c 
o 

u 

4J 

c 

<D 
O 

c 
o 
u 


200eV 


10 


-3 


■10 


0 
z(cm) 


10 


Figure  3.12 


Concentration  profiles  (p  =  0)  of  the 
molecular  hydrogen  distribution  obtained  from 
the  steady  state  solution  of  the  diffusion 
kinetics  system.   The  solutions  are 
normalized  by  Eq .  3.5. 
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Figure  3,13 


Concentration  profiles  (p  =  0)  of  the 
molecular  oxygen  distribution  obtained  from 
the  steady  state  solution  of  the  diffusion 
kinetics  system.   The  solutions  are 
normalized  by  Eq.  3.5. 
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Figure  3.14 


Concentration  profiles  (p  =  0)  of  the 

H  O  distribution  obtained  from  the  steady 

slale    solution  of  the  diffusion  kinetics 

system.   The  solutions  are  normalized  by  Eq 

3.5. 
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Compu  terAdaptation 

The  solution   of  the  system   of  equations   described  in 
chapter  III  is   based  on  a  finite   element  approach.    The 
basic  program  consist  of  the   IMSL  program  TWODEPEP  (1981). 
Extensive  input  and   output  routines  have  been   written  to 
apply  this  program  to  the  current  problem.    A  flow  diagram 
is  presented  in   Fig.   3.15.    The  program   is  listed  in 
appendix  B  and   the  main  input  routines  can   be  tailored  to 
the  time  dependent  or  steady  state  case. 

The  program  consists  of  five  main  input  subroutines. 

1.  COEFF   reads   in  scaling   parameters,    diffusion 
constants  and  miscellaneous  numbers. 

2.  REACIN  reads  in  the  species,  reaction  rate  and  yield 
parameters  and  sets  up  the  associated  matrices  which 
the  computational  subroutines  will  access. 

3.  MAIN  consist  of  a  number   of  general  eguations  which 
are  tailored  to  the  system  of  interest, 

4.  FS   is  the   main   computational  subroutine   which 
calculates  the   production  and  destruction   terms  in 
the  differential  system.   It  was  also  adapted  to  the 
kinetic  program  described  in  chapter  III. 

5.  FSD   calculates  the   partial   derivatives  of   all 
nonlinear  terms  with  respect  to  the  concentrations. 

Output  is  taken  care  of   by  subroutine  OUTPUT  which  produces 
a  grid   of  concentration   values  for   all  species   at  the 
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Figure  3.15   flow  diagram  for  diffusion  kinetic  code 
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requested  values  of  x  and  y.    This  subroutine  can  be  easily 
modified  to  process  the  iuformation   before  it  is  stored  on 
raaqnetic  tape. 

The  solution   algorithm  has  been  tested   extensively  by 
ir'iSL.    In  addition  the  procjram  was  tested  aqainst  a  simple 
one  diiuensional  system   similar  to  the  system   of  interest. 
Grid  dependent   tests  were   run  and  it   was  found   that  a 
minimum  of  130  finite  elements  was  needed. 


CHAPTEB  IV 
SESULTS 


JDiscussipn 
C  om£ u t a t i on al_S i m u la t i on 

The  results  of  the  calculations   can  be  divided  into  two 
broad  cateqories.    First,  the  calculation  of  initial  yields 
directly  resultiuq   from  the   deposition  of   the  incident 
electron  energy  in  the  medium.    Second,  the  calculation  of 
steady  state  concentrations  resulting   from  the  introduction 
of  electrons  into  the  medium  at  a  constant  rate.   The  system 
of  proqraais  which  are  used   to  calculate  these  quantities  is 
also  naturally   divided  into  tuese  two  categories.   These 
programs  are  listed  in  Fig,  4^1, 

The  programs   waich  are  used   to  calculate   the  initial 
yields  are  ELASTIC,  T0TH2C,  MESD,  GETMESD,  KC,  GETMC,  NLLSYS 
and  YIELD.    ELASTIC  and  T0TIi20   are  service  prograas  which 
are  used  to   generate  cross  sections  for   both  elastic  aiid 
inelastic  collisions.   In  addition,  they  also  produce  graphs 
ox  these  cross  sections.    MESD  and   GETHESD  are  a  pair  of 
programs   which  follow   the  path   of   an  electron   which 
undergoes  only  elastic  scattering.     MESD  generates  a  data 
base  of  the  collisional  history   of  the  electron  and  GETMESD 
sorts   this   information    and   calculates   cuiaulative 
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Figure  4.1   Overview  of  computation 
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distributions  for  R,  Z  and  X  as  a  function  of  pathlenqtli  and 
incident  electron  energy.    These  distributions  are  fit  with 
analytic  functions   which  are   then  inserted   in  the   MC 
proqram.    The   MC  proqram  is   the  maior  proqram   in  this 
seqment  of  the  calculation.    This  program  tracks  electrons 
on  a  collision   by  collision  basis  as   they  undergo  elastic 
and  inelastic   scattering.    When  the  electrons   reach  the 
point  at   which  the  total   elastic  cross  section   is  much 
greater  than  the  total  inelastic   cross  section  the  elastic 
distributions  produced  by   MESD  are  used  to   calculate  the 
movement  of  the  electron  between  inelastic  collisions.    In 
addition,   when  the  electron  energy  is  less  than  the  lowest 
excitation  state  being  considered,   these  distributions  are 
used  to   follow  the   electron  as   it  approaches   thermal 
eguilibriuia.    GETMC   uses  tiie  data   generated  by   MC  to 
calculate   two   and   four   dimensional   yield   spectra, 
collisional  densities   and  theriaal   electron  distributions. 
The  NLL3YS   program  is  used  to   fit  the  yield   spectra  and 
thermal  electron  distribution.    It  consists  of  a  nonlinear 
leasx  square  fitting   algorithm  which  has  been   modified  by 
the  incorporation   of  a   a  normalization   criterion  which 
reguires  that  the  integral  of  the  function  being  used  to  fit 
the  data  over  all  space,  integrate  to  unity.    The  laultiple 
dimensional  integrals  are  numerically   computed  in  spherical 
coordinates  using   Gauss-Laguerre  guadrature  in   the  radial 
coordinate  and   a  simple  Newton-Cotes  composite   formula  in 
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tiie  theta  coordinate.  Finally,  the  YIELD  program  uses  the 
analytic  fits  generated  by  NLLSYS  and  the  cross  section  data 
from  TOTH20  to  calculate  yields  as  described  by  Eg.  2.14.  A 
variation  of  this  program  generates  spatial  data  sets  whicu 
are  used  to  produce  isoyield  and  isodose  curves.  The  yields 
are  then  used  in  the  diffusion  kinetics  calculation. 

The  programs  which  are  used   to  calculate  the  subseguent 
evolution  of  the  system  are  KINETIC,   DIFKT,   and  CONTOUR. 
KINETIC  is  a  program  which   calculates  at  the  time  evolation 
of  a  set  of  first   order  nonlinear  reaction  rate  eguations. 
This  program  is   used  to  study  various   reaction  mechanisms 
without  diffusion.    DIFKT   is  the  maior  program   in  this 
calculation.    As  input,   it  reguires  the  results  from  the 
YIELD  program,  the  reaction  mechanism  studied  in  the  KINETIC 
program  and   a  set  of   diffusion  coefficients  in   order  to 
solve   a   system   of  second   order   nonlinear   partial 
differential  eguations  which  describes  the  physical  chemical 
system  set   up  by  a  steady   beam  of  electrons   incident  in 
gaseous  ii^C, 

i52^te_Carlg_and_Yieid3 

The  solutions  to  tue   multiple  elastic  scattering  problem 
presented  in  chapter   II  were  obtained  for   electrons  with 
energies  from  .2  eV  to  50  eV.    It  is  of  interest  to  analyze 
theiie  solutions   in  terms  of   the  single   scattering  cross 
section  at  these  energies  and   look  at  the  energy  dependence 
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of   the  parameters   used   in   representinq  the   multiple 
scattering  distributions  (Fiqs.  4,2-4.7).    First,  we  should 
note   that   at   all  energies   HO   has  a   significant 
backscatterinq  component-     For  example,   at  50   eV  the 
scattering  cross   section  is   essentially  split   into  two 
components  with  no   scattering  at  90  degrees.     At  lower 
energies  (less  than  1  eV)  we  see  (Fig.   2.1)  that  the  single 
scattering  distribution  is  monotonically   decreasing  and  the 
relative   difference  in   magnitude   between  the   forward 
scattering   peak,  and   the   backward   scattering  peak   is 
approximately  a  factor  of  ten.   At  higher  energies  {15  eV  to 
50  eV)    the  probability  for   small  angle   scattering  has 
increased  as  expected,   bat  the  ratio  of  the  areas  under  the 
curve  from   90  to   180  degrees   and  0   to  90   degrees  is 
increasing.   Looking  at  the  cumulative  radial  distributions, 
F(n)  in  Figs-   2.6-2.8,   we  can  see  that  as  the  pathiength 
increases,  the  distribution  widths  increase.    Furtherfiiore, 
as  the   electron  energy  increases  the   distributions  spread 
out,   the  increased  dispersion  correlating  witn  the  increase 
of  the  forward  scattering  aiaximuffl-    As  the  energy  continues 
to  increase,   these  distribution  widths  decrease  because  the 
ratio  of   backscattering  tc  forward   scattering  increases. 
These  patterns  are  expressed  by   the  parameter  dependence  on 
energy  as  depicted  in  Figs.   4-2  and  4.3-    The  cumulative 
longitudinal  distributions   behave  similarly.     The  laost 
apparent  feature  of  the  polar   angle  distributions  is  their 
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relative  uniformity  for  pathlenqths   greater  than  four  raeau 

free  paths,  (Fiqs.  2.12-2.14)  at  all  energies.   This  feature 

simplifies  the  application  of  the   HESD  results  in  the  Monte 

Carlo  transport  code  since  for  pathlenqths  greater  than  this 

it  is  not  necessary  to  invert  this  function.    In  practice, 

the  MESD  results  are  not  used  above  an  energy  v/here  there  is 

a   20r3  probability   of  having   fewer   than  four   elastic 

collisions   between  inelastic   collisions.    The   energy 

dependence  of  the   polar  angle  distribution  is  graphed  in 

Fig.  4.4. 

The  results  of  the  Konte   Carlo  calculation  are  presented 

in  condensed  form  by  the  analytic  representations  of  the  two 

and  four   dimensional  yield   spectra.    Samples   of  these 

functions  are  presented  in  Fiqs.  2.15-2.17.    The  functional 

representation  of  this  information   originated  in  the  papers 

by  Green,   Garvey  and  Jackman  (1977)   and  Green  and  Sinqhal 

(1979)   for  the   two  and  four  dimensioiial   yield  spectra, 

respectively,   and  were  used  as  recently  as  1933  in  a  study 

of  electron  deposition  in  neon  by  Singhal,  Rio,   Schippnick 

and  Green  (1983).   In  this  study  the  representations  for  the 

four  dimensional  spectra  for  H  0   was  better  fitted  with  the 

2 

new  functional  forms  presented  in  Tables  2.10  and  2.11. 

Two  diaiensional  yield  spectrum  are  qraphed  in  Fig.   2.10. 
It  is  interesting  to  compare  the   spectra  for  HO  with  those 
presented  for  some  other  gases  in  a  paper  by  Green,   Jackman 
and  Garvoy  (1977).    Some  of   these  have  been  reproduced  in 
Fig.  4.8. 
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Two  dimensional  yield  spectrum  for  various 
atoms  and  molecules.   (Green,  et  al.,  1977) 
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By  and  larqe,    ail  the  spectra  show   the  same  overall 

behavior  except   for  a   prominent  difference   in  the   low 

spectral  enerqy  behavior.    The  difference  can  be  attributed 

to  the  structural  variation  of  the  scattering  conponent,  for 

example  the  difference  between  an   atom  and  molecule.    Por 

both  H  0  and   CO   the  density  of  vibrational   states  at  low 
2  2 

enerqy,   as  plotted  in  Fiq.2.2  for  H  0  qives  rise  to  an  area 

2 

of  increased  collisions  at  low  enerqy.    Diatomic  molecules 

show  less  of   this  behavior  in  their   yield  spectra  because 

they  have   only  one  low   enerqy  vibrational   state  whereas 

their   spectra  are   comparable  to   those   of  atoms   (see 

Fiq.4,8).    In  addition,  a  two  dimensional  spectrum  for  E   0 

2 
calculated  in  1977   is  shown.    The  low   enerqy  resonance 

behavior  for  the   vibrational  states  was  not   known  at  that 

time,  and  by  comparinq  these  spectra  (Fiq.4.3)   to  those  in 

Fiq.   2.15  it   is  obvious  that  their   contribution  to  the 

spectral  distribution   is  shown  by   the  rapid   increase  in 

collision  density  at  low  enerqy  (see  Fiq.  2.2). 

The  functional   forms  for   the  four   dimensional  yield 

spectra  represent   an  improvement  in  the   representation  of 

the   physical  attributes   of  the   deposition  pattern   of 

electrons  in  a  qaseous  medium.    The  functions  used  to  model 

the  two   and  four  dimensional   yield  spectra   retain  their 

overall   attributes  as   formulated  by   Green  and   Jinqhal 

(1979),    That  is,  they  are  separable  in  spectral  enerqy  and 

spatial  coordinates,   the  four  dimensional  yield  spectra  are 
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constrained  to  inteqrate  to  the  two  dimensional  spectra,  and 
the  spatial  dependence   is  assumed  to  be   distributed  over 
three  spectral  energy  regions.    However,  the  current  model 
has  replaced  the  spatial  functions   with  functions  that  are 
based  on  concentric   ellipses  (see  Tables  2.10   and  2.11). 
Unlik:e  the  previous  functions,  these  are  continuous  in  the  z 
coordinate.    This  shape  is   similar  to  that  experimentally 
observed  by  Grun  (1957).   In  addition  the  density  dependence 
associated  with  the  four  dimensional  yield  spectrum  has  been 
corrected.    In  the  previous   representation  the  integrated 
yields  were   density  dependent   but,   as   can  be   easily 
demonstrated  by  looking  at  a   proiection  of  a  electron  path 
on  a  surface  which  expands  as   one  over  the  density,   it  is 
seen   that  while   the  spatial   pattern   may  expand   with 
decreased  density  the  number  of  total  collisions  and  in  fact 
the  spatial  pattern  remains  the  same.   This  concept  has  been 
incorporated  into   the  spatial  functions.     Finally,   as 
previously  mentioned,  the  four  dimensional  yield  spectra  are 
constrained   to  inteqrate   to  the   two  dimensional   yield 
spectra.    This   is  accomplished   by  requiring   that  the 
constants  A  observe  the  normalization  relationship  expressed 
in  Eg.  2.11.    Unfortunately,  the  spatial  integrals  are  not 
analytic  and  the  integrations   are  performed  numerically  as 
described  in  Chapter  IV.    Furthermore,   since  the  spatial 
parameters  P,   *,   and  z^  are  dependent  on  incident  energy, 
the  normalization   constants  A  are   also  dependent   on  the 
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incident  energy.    These  constants  are  tabulated  in  Table 
2.13. 

One   of  the   aost  interesting   results   of  using   the 
concentric   elliptic   functions  to   model   the   spatial 
dependence  of  the  yield  and   deposition  spectra  is  that  the 
incident  energy  dependence  of  the   spatial  parameters  can  be 
related  to  an  analytic  form  similar  to  that   used  for  the 
total  cross  section.    That  is,  3  can  be  thought  of  as  being 
related  to   the  total  cross   section  since   the  collision 
density  would   fall  off  slower  for  a  particle   which  had 
collisions  spread   far  apart.     Similarly  the   z  offset 
increases  with   an  increase  in  the   collisional  pathlength 
since  it   is  proportional   to  one   over  the   total  cross 
section.    The  A  parameter  controls  the  eccentricity  of  the 
ellipses.    In  general  it  is  seen  to  behave  in  a  predictable 
manner. 

The  yields  for   most  of  the  states  listed   in  Tables  2. 4 
and  2-5  are  calculated  using  Eg,  2.14  and  presented  in  Figs. 
4.9-4.13-   The  G(#/eV)  values  are  seen  to  become  constant  at 
high  energy.    The  value  of   W ( eV/ip)   for  the  ionization 
states  versus  incident  electron  energy   is  presented  in  fig. 
4,14  and  compared  with  experimental  data.    The  variation  of 
this  guautity   with  incident  energy   has  been   verified  by 
numerous  experiments  over  the  years   and  is  accounted  for  by 
the  properties   of  the   yield  spectrum   {Green, Garvey  and 
Jackman ,  1 977) . 
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Figure  4.14 


v;  (#ip/eV)  vs.  incident  electron  energy. 
Experimental  results  from  Klots  (1968)  (  • 
and  Berger  and  Seltzer  (1964)   (A). 
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Looking   at  the   yield   contours   presented  in   fiqs. 

2.18-2.23  it   is  easy   to  see  that   the  pattern   for  each 

species   presented   is  characteristic   of   the   physical 

parameters  of   the  collisional  cross  section.     For  this 

reason  we  see  that  OH,   which  can  be  produced  by  states  with 

lover  energy  thresholds  than  H  0+   can,   has  a  more  diffuse 

2 

pattern  than  H  0+.     This  is  because  the   electron  energy 
2 

diminishes  as  it  moves  away  from  the  incident  point.    Also, 
as   the  incident   electron  energy   increases  the   thermal 
electron  distribution   becomes  more  spherical.     This  is 
because  the  electrons,   after  as   few  as  four  collisions  at 
low  energy,   lose  the  knowledge  of  their  initial  orientation 
(see  Figs.  2.12-2.1^4). 

Pif f usion_Kinetics 

The  results  from  the  KINETIC   and  DIFKT  programs  provide 

us  witii  a  unigue  view   into  the  physical-chemical  mechanisms 

involved  in  the  picosecond  region.   Refering  to  Fig.   3. 1  we 

can  see  that   the  ions  which  are  initially   produced  by  the 

incident  electrons   react  very  quickly  and   are  neutralized 

essentially  within   a  microsecond.     The  H  O   clustering 

2 

reactions   enable  the   proton  to   remain   stable  for   a 
relatively  long  time.    This  is  because  of  the  shielding  it 
receives  from   the  polar   molecules.    The   electrons  are 
finally   seen  to   combine  with   these  positively   charged 
species  until  only  the  incident   electrons  remain.    At  the 
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same   time   these   species   are   beiaq   depleted   the 
concentrations  of  the  radical  H,  0,   and  OH  are  increasinq. 
The  chemistry  which  involves  the   OH  radical  is  completed  in 
one  second.    In  approximately  two  seconds  the  H 2O  chemistry 
has  essentially  stabilized  and  we  are  left  with  H 2#  0 2,    H 2O 2 
and  a  small  amount  of  HO2.    However,   we  should  note  (Fiq. 
3,1)  that  hydrogen  has  not  yet  reached  steady  state. 

DK   calculations  were   performed   for  electrons   with 
incident  enerqies  of  200  and   1000  eV.    These  tvo  enerqies 
are  representative  of  the  types  of  deposition  patterns  which 
can  be  seen  in  the  energy  range  from  100  eV  to  10  keV.    The 
200  eV  distribution  patterns  for   all  the  primary  yields  are 
much  tighter  than  those  for  1  keV  (see  Fiqs,  2,18-2,23).   As 
will  be  seen,    this  creates  a  dramatic   difference  in  the 
steady  state  concentration  distributions. 

We   beqin  the   discussion   of   these  calculations   by 
analyzinq  the  results  for  the  1  keV  electron.    The  primary 
species  involved  can  be  divided  into  three  major  types.   The 
electrons,  the  posixive  ions  and  the  neutrals  OH,  H,   C  and 
H  9.    As  mentioned  previously,   the  electrons  are  the  only 
species  considered  which  are  not   strictly  conserved  in  the 
system.   In  the  initial  formulation  of  the  problem  they  were 
introduced  into   the  system   by  two   methods  in   order  to 
compensate  for  this  imbalance.    However,  it  was  discovered 
that  neither  fix  mentioned   was  necessary.    The  resultant 
steady  state  distribution  for  the   electrons  is  presented  in 
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Fiq.  3.2.    The  distribution  is  presented  in  profile,  as  all 
results  are,  and  shows  that,   as  expected,   near  the  oriqin 
the  steady  state   distribution  is  dominated  hy   the  yields. 
This  region  extends  to   approximately  seven  centimeters  for 
the  density  at   which  the  calculation  was  made   and  can  be 
compared  directly  to   the  yield  profile  presented   in  Fiq, 
3.3,         In  spite  of  the   fast  reaction  rates  associated  with 
the  electron   their  numerical  abundance  prevents   them  from 
beinq  significantly  depleted.    Thus,  those  electrons  which 
do  not   react  with  a   corresponding  ion  quickly   beqin  to 
diffuse   in  the   classical  diffusion   liait.    This   is 
demonstrated  by  the  1/r  shape   of  the  profile  for  distances 
greater  then  seven  centimeters-     Furthermore,   we  see  by 
lookinq  at  the  profiles  alonq  both  the  p   and  z  axes  that  the 
distribution   is  almost   spherical  except   for  a   slight 
asya^raetric  component   in  the  forward  z   direction.    This 
asymmetry  (Fig.  3^2)         is  due  to  the   preference  of  the 
incidant  electrons  to  scatter  in  the  forward  direction. 

The  H2O+  ion  distribution  (Fiq.    3.4)   drops  off  very 
quickly.    This   is  because   the  maior   pathway  for   the 
destruction  of  this   ion  is  its  reaction  with   H2O  and  the 
reaction  rate  for  this  reaction  is  high.   We  see  in  Fig.  3.5 
that,  in  fact,   the  ion  reaction  is  so  guick  that  it  barely 
diffuses  before  reacting.   The  H2O+  distribution,  unlike  the 
electron  distritjution,  is  fairly  symmetric.   The  reaction  of 
the  H2O+  ion  leads  to  the  production  of  the  H3O+  ion.   Since 
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the  reaction  which  creates  this  ion   is  linear,   due  to  the 
constancy  of  the  H 2O  backqround,   the  distribution  for  this 
ion  is  almost  exactly  like  that  of  H 2O+ .    This  can  also  be 
attributed  to   the  fast  reaction   rate  of  H3O+   with  H2O, 
similar  to  its  precusser  fl  2O+ -    "^^^    remaininq  reactions  in 
this  subset  of  ionic  reactions,   as  presented  in  Table  3.2, 
continue  the  same  pattern.   That  is,  the  ions  form  complexes 
with   the  water   molecules  and   are   neutralized  by   the 
electrons  (Fiq.  3.6).    The  maior  effect  of  these  hydration 
reactions  is  to   delay  the  production  of   hydroqen  and  thus 
shift  the  hydroqen  concentration  away  from  the  oriqin.    In 
addition,   these  reactions  lead   to  decreased  production  of 
atomic  hydroqen  when   compared  with  the  direct   reaction  of 
the  U  2O+   i°i^  with   the  electron.    This   is  due   to  the 
exothermicity  of  the  charqe  neutralization  reactions  for  the 
hiqher  order   hydrates  which   in  turn   leaves  insufficient 
enerqy  to  dissociate  a  water   molecule  as  compared  with  the 
neutralization  reaction   of  the   hydroniuia  ion   (Anderson, 

1958)  . 

Next,  let  us  look,  at  the  steady  state  distribution  of  the 
hydroqen  atom.    Unlike  the   H  2O+  ^°^    previously  discussed, 
the  maior  source   of  these  atoms  is  not  solely   due  to  the 
dissociation  of  the  water  molecule   by  the  incident  electron 
but  also  due  to  the   neutralization  process.    We  therefore 
see  a  rather  extended  reqion  of  constant  concentration  level 
for  this  species  (see  Fiq.    3.7)   even  thouqh  the  primary 
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yield  falls   off  much  quicker   (Fiq,   3.8) .  The  large 

diffusion  coefficient   of  this  species  also   contributes  to 
this  behavior.    The  hydrcqen  molecule  distribution  follows 
the  same   pattern  as  that   of  atomic  hydroqen.     This  is 
because  the  major  pathway  for   production  of  this  species  is 
via  the  combination  of  atomic  hydroqen. 

The  OH   radical  exhibits   a  distribution   pattern  which 
shares  some  attributes   with  both  the  H2O+  ion   and  0  atom 
distribution.    For   instance,   the  OH   distribution  (see 
Fiq.  3-9)  beqins  to  decrease  almost  immediately  because  OH  is 
fairly  reactive.    However  we  can  see  that  the  distribution 
does  not   fall  off   as  rapidly   as  that   for  the   ions. 
Furthermore,   by  coraparinq  the  yield  profile  with  the  steady 
state  distribution  profile  (Fiq.   3.10)  we  can  see  that  the 
radical  moves  an  appreciable  distance  before  it  reacts.   The 
oxyqen  atom  distribution  (see  Fiq.   3.11)  is  similar  to  that 
of  the  OH  radical.    Essentially  its  fall  ofx  is  slower  than 
that  of  OH  because  Atomic  oxyqen   does  not  react  as  readily 
as  does  OH  and  taerefore  its   distribution  is  not  as  sharply 
peaked.    The   molecular  oxyqen  steady   state  distribution 
looks  similar  to  that  of   molecular  hydroqen  except  that  the 
plateau  reqion  is  not  as  flat.    This  is  expected,   to  some 
extent,   since  O2  and  H2  occupy  somewhat  of  the  same  role  in 
the  cheaical  mechanism.    That  is,  they  are  botn  relatively 
stable  cnemical  species.   The  final  products  of  interest  are 
H2O2  ^^'^      ^^02*    ^^^^  hydroqen   peroxide  molecule   is  not 


120 
produced  in  any  appreciable  amount.    This  is  not  unexpected 
since  the  kinetic  model  gives  an   indication  tiat  it  is  not 
stable  enouqh  to  exist  for   any  significant  length  of  time. 
On  the  other  hand,   H2O2   appears  relatively  stable  in  the 
kinetic   calculation  and   we  see   that  its   concentration 
distribution  is   significant.    The   peak  value   of  this 
distribution  is  of  course  much   smaller  than  that  of  either 

H2  or  ©2. 

Eefering  to  Figs.   3.2,  3.4,   3.7  and  3.9  we  can  observe 
the  critical   differences  in  the  steady   state  distribution 
patterns  for   the  two  incident  electron   energies  studied- 
The  low  energy  incident  electrons  give   rise  to  a  much  more 
compact  region  of   positive  ions  which  in   turn  deplete  the 
thermal  electron  population.    This  inner  core  of  positive 
charge,   while  present   at  higner  energy,   does   not  peak 
dramatically  enough  to  deplete   the  thermal  electrons.    In 
general  this  pattern   of  high  peak  values  with   a  fast  fall 
off  is   seen  for  all  the   species  created  directly   by  the 
incident  low  energy  electrons.   The  distribution  patterns  of 
the  priaary   species  produced   by  the   1000  eV   incident 
electron  show   a  more  gradual   transition  region   to  pure 
diffusion.    While  there  is  a  significant  difference  in  the 
distribution  patterns  of   the  priaary  species  for   the  two 
cases  studied,    those  species  which  are   essentially  end 
products   of  the   reaction  mechanism   are   seen  to   have 
distributions  which  are  similar  in   shape  (see  Fig-   3.12, 
3.  13  and  3.14)  . 
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The  results  from  the  diffusion   kinetic  proqram  DIFKT  are 
presented  in  Fiqs.   3.2-3-14  for   200  and  1000  eV  electrons 
incident  at  the  rate  of  10(12)   per  second.    Cotnparinq  the 
profiles  of   the  steady  state   concentration  distributions 
with  the  yield  profiles,   we   see  that  the  distributions  of 
those  species   directly  produced  by  the   incident  electron 
exhibit  an  appreciable  relaxation   of  the  system.    Whereas 
the  yields  fall  off  dramatically  with  distance,   the  steady 
state  concentrations   fall  off   much  more   slowly.    The 
electron,   whose  diffusion  coefficient  has  been  estimated  at 
approximately  five   times  that  of   H  or  H2 ,    travels  the 
farthest  distance   before  it   beqins  to   be  depleted   by 
interactions  with  the  positive   ions.    Molecular  hydroqen, 
whicn   is  one   of  the   major  products   of  tae   reaction 
mechanism,   is  seen  to  remain  constant  for  a  relatively  lonq 
distance.    An  analoqous  behaviour  is  also  exhibited  in  the 
kinetic  results   presented  in  fiq.    3,1.    This   can  he 
attributed  to   the  relative   stability  of   tnis  molecule. 
Furthermore,  we  note  (Fiq.   3.1)  that  li,  is  the  last  species 
to  stabilize.     Therefore  wo  would  expect   the  molecular 
hydroqen  distribution   profile  to  extend  further   out  than 
that  of  some  of  the  other  species,   like  H^  0+  or  OH,   which 
react  very  quickly. 

The  picture  of  the  diffusion   kinetic  system  that  we  see 
is  essentially  one  of  translatinq  the  time  coordinate  in  the 
time  dependent  kinetic  simulation   (Fiq-   3-  1)   into  space 
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where   we  have   an   extended   source.    The   resultant 
concentration  profiles  are  essentially  spherical   and  fail 
off  gradually  in   comparison  to  the  source   term  (yields). 
The  electrous  form  an  outer  layer  of  charqe  around  the  ions, 
which  of  course  remain  nearer  the  oriqin,   because  of  their 
slower  diffusion  rate.   The  radicals,  for  example  OH,  have  a 
profile   which  extends  further   out.    Finally,    the 
distribution  profiles  for  H  ,   0^  and  fi2°2 '  ^''^^^   '^^^   ^^'® 


m 


ost  stable  species,  extend  out  the  farthest. 


Limitations 
Let   us  discuss   the  iimitatious   of  the   raatheaatical 
techniques  used  in   this  study  to  describe   the  effects  of 
irradiating  H  0  vapor  with  electrons.   Foremost  is  the  yield 
spectrum  which   characterizes  the  immediate  iiapact   of  the 
electron  on   the  medium.    The  first   maior  consideration, 
especially  in  a  gaseous  medium,   is  that  the  density  of  the 
medium  be  sufficiently  high  that   the  radiation  not  diainish 
it  to  an  appreciable  extent.    This  is  necessary  because  the 
yield  spectrum  is   dependent  on  the  physical   properties  of 
the  medium.    Since  we  have  chosen  to  study  the  steady  state 
case  for  the  diffusion  kinetic   system  it  is  also  necessary 
to  consider  when   this  situation  might  break   down.    The 
system  looked   at  assumes  a   medium  which   is  essentially 
infinite  compared  to  the   spatial  distribution  of  products. 
Thus  we  have   not  taken  into  account   the  finite  boundari^is 
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which  would  be  associated  with   a  containment  vessel  for  the 
medium  which  of  course  could  be  bioloqical.    Furthermore, 
the  chemistry  associated   with  the  system  (see   Fiq,   3*1) 
implies  that   the  ions  are  depleted   in  the  time   frame  of 
10  (-6)   seconds.    Therefore  the  rate  of  input  of  electrons 
must   be  greater   than  one   over  this   number  to   avoid 
calculations  with  less  than  one  electron*    This  may  happen 
since   we  are   using   a   continuum  description  in   the 
mathematical  formulation  of  the   diffusion  kinetics  physical 
chemistry. 


CHAPTEE  V 
CONCLUSIONS 


In  the  early   fifties  a  qroup  of  scientists   set  in  Hiqh 
Park, Illinois  to  discuss  the  physical  mechanisms  involved  in 
the  interaction  of  radiation  with   a  bioloqical  cell,    la 
1964  and  aqain  in  1974  in  Airlie,  Virginia  the  problems  were 
addressed  by  an  invited  qroup   of  scientists  froK  the  fields 
of  bioloqy,  chemistry  and  physics.   Durinq  the  1974  meetinq, 
the  concept  of  following   the  radiation  deposition  mechanisa 
from  initial   principles  was   seriously  discussed   for  the 
first  time.    One  of  the   founders  of  this  field,   Robert 
Platzman  (1974),   said   in  a  review  article   oa  early  time 
chemistry  froia  this  meetinq  that  "at   a  meetinq  of  this  type 
20  years  aqo  and  at  another  10  years  aqo,  I  encountered  some 
of  the   same  problems  I  have   touched  on  here,   and   I  am 
willinq  to  bet  anyone  that  if  I  survive  another  10  years,   I 
will  find  that   some  of  the  complications   I  have  stressed 
have  survived  also."   In  the  last   ten  years  the  research  in 
this  area  has  literally  exploded.     The  latest  meetinq  of 
participants  in   this  interdisciplinary  field  was   held  at 
Gettysburq,   Pennsylvania  in  1932  in  what  was  the  fifteenth 
annual   meetinq  on   "The  Bioloqical   Effects  of   lonizinq 
Radiation-   Radioloqical  Physics".    This  meetinq  led  to  a 
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special   "Workshop  on   the   Interface  Between   Radiation 
Chemistry  and  Radiation  Physics"  held  later  iu  the  same  year 
at  Arqonne  National  Laboratory. 

As  mentioned  in  the  introductory  chapter,   the  questions 
raised  thirty  years  aqo  are  finally  beinq  answered.   In  this 
dissertation   we  have   brouqht   toqether  information   and 
techniques  from   many  fields  in  order   to  qive  one   of  the 
first  complete  descriptions  of  the   physical  and  early  time 
physical-chemical  raechanisms   involved  in  the   radiation  of 
bioloqical  systems.   Thus  we  have  followed  the  deposition  of 
enerqy  by  an  electron  as  it  produced  initial  yields  and  have 
tracked  these  species  as  they  reacted  and  diffused.   The  sat 
of  proqraas  written  to  describe  this  system  provides  us  with 
an  apparatus  to  study  various   physical  and  caemicai  effects 
in  a  very  detailed  manner.     For  example,   the  effect  of 
varyinq  the  dose  rate  or  enerqy   of  an  electron  beam  can  be 
seen  from  a  new  perspective.    It  is  also  possible  to  study 
the  effects  of  addinq  trace   amounts  of  various  chemicals  to 
the  system  and   lookinq  for  a  unique  mechanism  which  would 
modify  the  effects  of  the  radiation. 

As  a  byproduct   of  these  calculations  we   have  computed 
yields  for  the  various  individual   states  and  seen  that  they 
are  similar  to  those  observed   in  experiments.    It  follows 
that  the  ev/ion  pair  calculated  is   also  quite  close  to  tae 
observed  value.    Furthermore,   we   are  able  to  construct 
isoyields  or,   of  more  importance  to  the  practicinq  medical 
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physicist,   isodose   curves  for   a  moaoeiierqetic   beam  of 
electrons.    It  is  a  straightforward  procodure  to  qeneralize 
these  computations  for  a   spectral  and  spatially  distributed 
beaiR  of  electrons. 


APPENDIX  A 
MONTE  CARLO  PBOGRAtJ 


The  proqraa   used  to  qenerate   the  path  of   an  electron 

inelastically  scattering  is  MC.   This  program  was  originally 

written  by  Brinkman  in  1970  and  essentially   rewritten  by 

Jackman  (1978).     The  program   has  since   bean  slightly 

modified  by  Singhal  (Sinqhal  et  ai.   1980,   Sinqhal  et  al . 

1983).    la  the  course  of  calculating  the  yield  spectrum  for 

H  0  this  program  was  again   revised.    Approxiuiately  70,<  of 

2 
the  program  was  rewritten;  30%   of  this  was  new,  reflecting 

changes  in   both  the   logic  of  the   program  in   order  to 

decrease  its   computational  time   and  to   accommodate  the 

necessity  of   calculating  the   collisional  history   of  the 

electrons  as  they  approach  thermal  equilibrium.    There  were 

some  problems  associated   with  this  program  prior   to  these 

modifications.   During  the  research  on  Neon  {Singhal  et  ai., 

1983)   the   probabilities  for   inelastic  collisions   were 

incorrectly  scaled.     The  allowed   and  forbiddfn   state 

inelastic  scattering   probabilities  were  scaled  only   by  a 

proportionality  factor  above  200  eV,   whereas  the  forbidden 

probabilities  fall  off   much  faster  than  the   allowed  with 

increasing  energy.    Additional  problem   areas  included  the 

omission   of   scattering  for   the   electron   immediately 
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followinq  the  MESD  collision,    cross  correlations  between 
MESD  parameters  because  the  same   random  number  was  used  for 
generation  of  uncorrelated  variables  and  finally  the  ratio 
of  elastic  to   inelastic  cross  sections   used   in  the  MESD 
subroutine  was  incorrectly  calculated   since  the   elastic 
cross  section   changes  for  large  pathlengths   (see  chapter 
II) .    In   order  to  clarify   this  problem  with   the  MESD 
subroutine  in  tlie  in  the   Honte  Carlo  transport  calculation 
it  is  helpful  to  refer  to  Fig.  K.2.      Originally  when  an  MESD 
calculation  was   applied  between   the  origin   (relatively 
speaking)   and  the  next  inelastic  collision  (the  second  set 
of  axes)    the  inelastic   scattering  was   ignored.    The 
inelastic  scattering  algothrim  introduced   into  the  current 
program  is   diagramed  in  Fig.    A.  3.    The   progra^a  was 
completely  annotated   and  the  input   was  modified   in  the 
interest  of  siiaplif ication.    In  addition,  a  logical  pathway 
was  created  which  allowed  us  to  follow  the  electrons  as  they 
elastically   scatter  iselow   the   lowest  inelastic   state 
considered.    The   electrons  were  then  tracJced   until  they 
approached  thermal  eguiliJarium, 

Following  these   changes,   the  program  was   checked  by 
running   a  few   test  electrons   through  and   calculating 
different  parts  of   the  path  with  built   in  diagnostics  aui 
coiaparing  them  with  hand  calculations.    The  diagram  in  Fig. 
A.I  shows   the  new  logic  flow   for  the  l\Z      proqram.    The 
program  was  found   to  be  approximately  50  ;i   faster  than  the 
older  version. 
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The  companion  proqrain  to  MC,    GETMC  was  also  modified. 
The  changes  reflect  those  made  in  the  HC  program  in  order  to 
follow   the  low   energy  electrons.     However,   a   few 
modifications  were  made  in  both  the   input  and  output  of  the 
program.    The  new  output  requires  some  discussion,  since  it 
enabled  a  longtime   embedded  problem  with  this   code  to  be 
found.    The  GETMC  code  essentially   calculates  the  two  and 
four  dimenional  yield   spectruas.    In  doing  this   it  sums 
collisional  events  in  energy  and  spatial  bins  and  normalizes 
the  results.    Unfortunately,  the  information  describing  the 
number   of  collisions   which  occur   in  each   bin  is   not 
displayed.     In  reworking   this  program   it  was   deemed 
necessary  to  output   this  information  in  order   to  fit  the 
yield  spectrum  data   only  where  the  number   of  collisional 
events  was   statistically  meaningful.    This   uncovered  a 
computer  error  which  caused  bins   with  no  collisional  events 
to  have  the  number  of   collisional  events  displayed  from  the 
previous  calculation.    Finally,   this  code  was  modified  to 
calculate  the  ion  pairs  per   eV  directly  from  the  tabulated 
data  instead  of  using  the   yield  calculation  as  discussed  in 
chapter  IV.    In  closing  this  appendix,   it  should  be  made 
perfectly  clear  that  due  to  the  nature  of  the  yield  spactrua 
most  of  the   errors  mentioned  will  not  have  a  significant 
effect  on  the  prior  results. 
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Figure  A.l   Flow  diagram  for  Monte  Carlo  transport  code. 
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Figure  A. 2   Coordinates  for  multiple  elastic  scattering 
calculation. 
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Figure  A. 3   Ionization  scattering  algorithms  E  is 
primary  electron  energy  and  T  is 
secondary  electron  energy. 


APPENDIX  B 
DIFFUSION  KINETIC  PROGEAl 


This  appendix  consist  of  the  source  code  which  interfaces 
with  TWODEPEP,   in  the  program   DIFKT  as  presented  in  Fiq, 
3.6. 
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